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The conditions for both the stability and the breakdown of the topological classification of gapped 
ground states of noninteracting fermions, the tenfold way, in the presence of quartic fermion-fermion 
interactions are given for any dimension of space. This is achieved by encoding the effects of 
interactions on the boundary gapless modes in terms of boundary dynamical masses. Breakdown 
of the noninteracting topological classification occurs when the quantum nonlinear sigma models 
for the boundary dynamical masses favor quantum disordered phases. For the tenfold way, we 
find that (i) the noninteracting topological classification Z 2 is always stable, (ii) the noninteracting 
topological classification Z in even dimensions is always stable, (iil) the noninteracting topological 
classification Z in odd dimensions is unstable and reduces to Zjy that can be identified explicitly for 
any dimension and any defining symmetries. We also apply our method to the three-dimensional 
topological crystalline insulator SnTe from the symmetry class AII-l-77, for which we establish the 
reduction Z —Zs of the noninteracting topological classification. 


I. INTRODUCTION 

Topological insulators (TIs) and topological supercon¬ 
ductors (TSs) of noninteracting fermions are character¬ 
ized by topological numbers (Z or Z2) that encode the 
non-trivial topology of the occupied single-particle wave 
functions and are accompanied by gapless excitations 
that are localized along any boundaryfii^ The integer 
quantum Hall effect (IQHE) is characterized by the Hall 
conductivity quantized by the integer v = 1,2 ,-•• in 
units of e^/h. The topological integer v counts the num¬ 
ber of extended chiral edge modes propagating at the 
boundary of the sample. The Z2 topological insulator is 
characterized by the parity of the number of Kramers’ 
doublets of extended boundary modes. Together with 
polyacetylene and a two-dimensional p -\- ip supercon- 
ductor^i^ both instances are now understood to be non¬ 
trivial entries in the periodic table (i.e., the tenfold way) 
for noninteracting topological insulators and supercon¬ 
ductors i^“— 

The gapless modes appearing at the boundary in the 
IQHE are robust to both elastic and inelastic scattering 
resulting from one-body impurity potentials and many- 
body electron-electron interactionsSimilarly, the gap¬ 
less modes in the Z2 TIs are immune to both backscat- 
tering resulting from one-body impurity potentials and 
many-body electron-electron interactions, provided time- 
reversal symmetry (TRS) is neither explicitly nor spon¬ 
taneously brokenfiSrJ^ 

Given the robustness to many-body fermion-fermion 
interactions of the edge states in the IQHE, it was a re¬ 
markable observation made by Eidkowski and Kitaev in 
2010 that it is is possible to gap out eight Majorana zero 
modes localized at the end of a one-dimensional topo¬ 
logical superconducting wire through many-body inter¬ 
actions without closing the spectral gap in the bulk i^'^d^ 


In the terminology of the tenfold way^^— it was demon¬ 
strated in Refs. [Tj and [H that the Z topological classi¬ 
fication for the noninteracting one-dimensional symme¬ 
try class BDI, when interpreted as a superconductor, is 
(i) unstable to quartic contact interactions that neither 
break explicitly nor spontaneously the TRS, and (ii) this 
instability reduces the noninteracting topological classi¬ 
fication Z to Zg. 

Subsequently, noninteracting two-dimensional topolog¬ 
ical crystalline superconductors (TCSs) from the sym¬ 
metry class DHI-|-i? (where “J-R” indicates the pres¬ 
ence of an additional reflection symmetry) and three- 
dimensional topological superconductors from the sym¬ 
metry class Dili were shown in Refs, [l^ and [13 and 
Refs, [isl - l^ to display the reduction patterns Z —>• Zg 
and Z —)• Zj^g, respectively, when perturbed by quar¬ 
tic contact interactions that neither break explicitly nor 
spontaneously the defining symmetriesi^ The reductions 
Z —>■ Z4 and Z —> Zg for the three-dimensional symmetry 
classes Cl and AHI were obtained in Ref. [IJ. 

We present in Sec. m a method that allows to de¬ 
rive the reduction pattern of all noninteracting topologi¬ 
cal insulators and superconductors without and with re¬ 
flection symmetries for any dimensionality d of space in 
the presence of quartic contact interactions that neither 
break explicitly nor spontaneously the defining symme¬ 
tries. This method relies on the topology of the classify¬ 
ing spaces from K-theory. It extends the applicability of 
K-theory for obtaining the tenfold way of noninteracting 
fermions^i^ to obtaining the breakdown of the tenfold 
way induced by interactions. 

This method is applied first to the breakdown of the 
tenfold way in Sec. IHIl ^ In doing so, we prove the fol¬ 
lowing properties that we report in Table U 

1. All Z2 entries of the periodic table irrespectively 
of the dimensionality of space are stable to quartic 
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contact interactions. 

2. All Z entries of the periodic table when the dimen¬ 
sionality of space is even are stable to quartic con¬ 
tact interactions. 

3. Only the Z entries of the periodic table when the di¬ 
mensionality of space is odd are unstable to quartic 
contact interactions with a reduction pattern that 
is computed explicitly and shown to break the Bott 
periodicity of two for the complex symmetry classes 
and of eight for the real symmetry classes. 

This method is then applied to the three-dimensional 
topological crystalline insulators (TCIs) from the sym¬ 
metry class AII-|-i?, which are of relevance to SnTe, in 
Sec.IVl We show the reduction Z —Zg in the presence 
of quartic local fermion-fermion interactions. 

The strategy that we use to study the robustness of v 
boundary modes to quartic contact fermion-fermion in¬ 
teractions is inspired by the (unpublished) approach pio- 
neered by Kitaev in Refs. [IM andlH, see also Ref. IH. It 
consists of three steps. First, a noninteracting topologi¬ 
cal phase is represented by the many-body ground state 
of a massive Dirac Hamiltonian with a matrix dimen¬ 
sion that depends on v. Second, a Hubbard-Stratonovich 
transformation is used to trade a generic quartic contact 
interaction in favor of dynamical Dirac mass-like bilin- 
ears coupled to their conjugate fields (that will be called 
Dirac masses). These dynamical Dirac masses may vio¬ 
late any symmetry constraint other than the particle-hole 
symmetry (PHS).— Third, the v boundary modes that 
are coupled with a suitably chosen subset of dynamical 
masses are integrated over. The resulting dynamical the¬ 
ory on the (d — l)-dimensional boundary is a bosonic 
one, a quantum nonlinear sigma model (QNLSM) in 
[{d — 1) -|- l]-dimensional space and time with a target 
space that depends on v. The reduction pattern is then 
obtained by identifying the smallest value of v for which 
this QNLSM cannot be augmented by a topological term. 
The presence or absence of topological terms in the rele¬ 
vant QNLSM is determined by the topology of the spaces 
of boundary dynamical Dirac masses, i.e., the topology of 
classifying spaces. Now, K-theory provides a systematic 
way to study the topology of the classifying space. Hence, 
this is why the same approach that was used to obtain the 
tenfold way of noninteracting fermions can be relied on to 
deduce a classification of topological short-range entan¬ 
gled (SRE) phases [also known as symmetry-protected 
topological (SPT) phases] for interacting fermions<^ 

Other topological phases are also interesting on their 
own right. For example, bosonic SPT (SRE) phases show 
many novel topological phases driven W strong interac¬ 
tions. They have been reviewed in Ref. 1^. The classifica¬ 
tion of bosonic SPT (SRE) states has been obtained by 
diverse approaches that include group cohomology^Si^ 
the K-matrix approach;^ enumerating surface topologi¬ 
cal order^ wire constructionsand so on. Topolog¬ 
ical order with long-range entanglement (LRE) is also a 


subject of intensive studies, which have relied on parton 
constructionstopological field theories— exactly 
soluble models)^^— and wire constructions— r— 

II. STRATEGY 

In this section, we present our strategy to obtain a 
topological classification for interacting fermions with 
gapped ground states as an application of K-theory to 
certain dynamical Dirac masses for boundary fermions. 

Noninteracting fermions always belong to one of the 
ten Altland-Zirnbauer (AZ) symmetry classes defined by 
the presence or absence of the following three symmetries, 
time-reversal symmetry (TRS), particle-hole symmetry 
(PHS), and chiral symmetry (CHS) (see Appendix lAl) . 
Within any one of these ten symmetry classes, the defin¬ 
ing topological attributes of noninteracting topological 
insulators and superconductors are shared by equivalence 
classes of Hamiltonians. Any two members within a topo¬ 
logical class can be deformed into each other by a smooth 
(adiabatic) deformation of the matrix elements of these 
Hamiltonians without closing the bulk energy gap. These 
equivalence classes are endowed with an Abelian group 
structure ©. Eor any given dimensionality d of space, 
topological invariants 0 are non-trivial for five out of the 
ten AZ classes. Specifically, three of the ten AZ classes 
support Abelian groups © = Z, while two of the ten AZ 
classes support Abelian groups © = Z 2 . The TRS, PHS, 
and CHS can be augmented by crystalline symmetries. 
Noninteracting fermions obeying crystalline symmetries 
can also be understood as realizing topologically distinct 
equivalence classes, i.e., topological crystalline insulators 
(TCIs),^ 

The topological classification with the Abelian group 
© for noninteracting TIs, TSs, or TCIs can break down 
in the presence of many-body interactions. Namely, an 
Abelian group 0,^^^ that encodes the topological equiv¬ 
alence classes of gapped ground states for interacting 
fermions can be smaller than © as a group (some quotient 
group of ©). 

In order to establish the instability of the noninter¬ 
acting classification of TIs, TSs, and TCIs, we choose a 
family of massive Dirac Hamiltonians, 

:= —iy^ 5 (g] 1 -I- m{x)f3 ig 11 (2.1) 

as representative single-particle Hamiltonians. Here, 
Dirac matrices a and P anticommute with each other 
and have the minimal dimension (rank) r^^^ under the 
symmetry constraints, i.e., is the minimal rank to 
realize a Dirac Hamiltonian of the form (12.11) . The di¬ 
mension of the unit matrix ]lisj^ = l,2,---. The integer 

G © is then related to the dimension r[v) = r^^^ v of 
the Dirac matrices that we choose. The question that 
we want to address is that of the stability or instabil¬ 
ity of the boundary states of a noninteracting TI, TS, or 


3 


TCI in the presence of many-body interactions that do 
not break the protecting symmetries of the noninteract¬ 
ing limit.— Here, whenever ly ^ 0, the extended single¬ 
particle boundary states are governed by the massless 
Dirac Hamiltonian 

o 

■==-i^ (2-2) 

which is obtained by introducing a domain wall in the 
mass m(x) along the a;^-direction that enters Hamilto¬ 
nian m- The Dirac matrices a (D 1 have a dimension 
r(i/)/2 that is half that of the bulk massive Dirac Hamil¬ 
tonian The dimension of the matrices a is r^^^/2. 

The breakdown (reduction) of the topological classifi¬ 
cation for noninteracting fermions takes place when the 
boundary states of the TIs, TSs, or TCIs can be gapped 
by many-body interactions that preserve their defining 
symmetries. By assumption, we consider many-body in¬ 
teractions that are weak relative to the bulk gap. If so, 
it is sufficient to treat the effects of many-body interac¬ 
tions for the massless Dirac fermions propagating on the 
(d— l)-dimensional boundary. To establish an instability 
of the noninteracting topological classification, we need 
not consider all possible many-body interactions. It suf¬ 
fices to establish that at least one family of strong (on the 
boundary) interactions implies the instability of the non¬ 
interacting classification 0 by gapping out all boundary 
Dirac fermions. To this end, we limit ourselves to contact 
interactions. 

Contact interactions are constructed from taking 
squares of local bilinears in the Dirac fermions. We have 
two options for these bilinears. The bilinear under con¬ 
sideration either commutes or anticommutes with the ki¬ 
netic contribution to the Dirac Hamiltonian. We shall 
call the latter option a Dirac mass. In this paper, we 
only consider the contact interactions obtained from tak¬ 
ing squares of those bilinears built out of Dirac mass ma¬ 
trices, for only these can gap the noninteracting massless 
boundary Dirac fermions in a mean-field approximation. 
Because we assume that the protecting symmetries forbid 
the presence of Dirac masses on the boundary that are 
consistent with the protecting symmetries, the only pos¬ 
sible Dirac masses induced by a mean-field treatment of 
a symmetry-preserving quartic interaction on the bound¬ 
ary must be odd under at least one of the protecting 
symmetries. We shall call such a boundary Dirac mass 
a boundary dynamical mass and label it with the Greek 
letter /?. 

We are thus led to consider the many-body interacting 
Dirac boundary Hamiltonian 


and 




(int)_ 


bd 


:= A 


1/3} • 


X 


l3"^(t,x) . (2.3c) 


We have chosen the real-valued coupling A with the di¬ 
mension of (length)"^ ^ to be independent of (3 for sim¬ 
plicity. This coupling constant is marginal in d = 2 and 
irrelevant when d > 2. (Of course, it can very well be that 
the set {/?} is empty. If so, we anticipate that 0 = 0;^^^ 
must hold. This is what happens for the strong topo¬ 
logical insulators in the symmetry classes A, D, and C 
when d = 2.) At this stage, it is convenient to treat the 
many-body Hamiltonian (12.31) with the help of the path 
integral 


Zy,^:= (2.4a) 

where the action in Euclidean time r is 

S^^:= I dr (2.4b) 

with the Lagrangian density 

/ltd := + • (2-4c) 

1 / 3 } 

The path integral is over Grassmann-valued Dirac 
spinors. 

We rewrite the quartic interaction terms by performing 
a Hubbard-Stratonovich transformation with respect to 
the bosonic fields 4>p conjugate to 4'f/3 4>, 

ZbdCxyp[vI/,vI/t,</.^]e-^^>T (2.5a) 

Here, the action in Euclidean time r is 

SU:= I dr Jd^-^xC'^^, (2.5b) 

with the Lagrangian density 

:= {d, + + (2-5c) 

m 

where we have introduced the dynamical one-body single¬ 
particle Hamiltonian 

^bd^"^(Da^) :=^m( 3:)+Xl2i/^'^/3(Da^), (2.5d) 

1 / 3 } 


(2.3a) 

where (the subscript “bd” stands for boundary) 

:= / d‘^-^x 'kt(f, x) T(t, x) (2.3b) 


under the assumption that the sign A > 0 corresponds 
to a repulsive interaction. In a saddle-point approxima¬ 
tion, the magnitude of the vector (p with the compo¬ 
nents (pp can be frozen both in imaginary time and in 
{d— l)-dimensional space. Fluctuations that change this 
frozen magnitude are suppressed by the second term on 
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the right-hand side of Eq. (I2.5cl) . We will restrict the 
set 1/3} to pairwise anticommuting Dirac mass matrices. 
If so, the direction in which the vector (p with the com¬ 
ponents (j)p freezes in the saddle-point approximation is 
arbitrary!^ Since fluctuations about this direction are 
soft, these are the Goldstone modes associated with the 
spontaneous breaking of a continuous symmetry. 

The effective low-energy theory governing the fluctua¬ 
tions of these Goldstone modes is obtained from a gradi¬ 
ent expansion of the fermion determinant 

Det + := 

It is captured by the partition function 

\d ~ ^(0^ - 1) e“‘^«NLSM-'S'top^ (2.7) 


after we have rescaled the vector cp so that it squares to 
one. The Euclidean action 

5'qnlsm = ^ (2.8) 

is the action of the quantum nonlinear sigma model 
(QNLSM) with the base space in space and 

time and the target space 

(29) 


with the integer N{v) counting the pairwise anticommut¬ 
ing Dirac masses that have been retained in the set {/3}. 
The effective coupling constant g is positive. The topo¬ 
logical term is present whenever any one of the ho- 
motopy groups 


(^7V(.)-l) ^ 




( 2 . 10 ) 


u)-i 


'd+l 


(^gN(v)-l 


is non-vanishingi^ (The reason why we ignore all topo¬ 
logical terms associated with non-vanishing homotopy 
group of order larger than d -|- 1 is that such topolog¬ 
ical terms would modify the local equations of motion 
derived from /S'qj^Lgjyj in a non-local way.) It signals the 
existence of zero modes of the Dirac Hamiltonian (I2.5dl) 
in the presence of topological defects in the order param¬ 
eter p. We expect that these zero modes prevent the 
gapping of the boundary Dirac fermions. We define the 
smallest value for the dimension v of the unit matrix 


1 in Eq. (12.21) for which 


TTo = 0 , 

TTi = 0 , 

TT, = 0 , 

TT^+l = 0. 


( 2 . 11 ) 


As all homotopy groups of the spheres are known, one 
may verify that 


d+l<iV(j.„i„)-l. (2.12) 

When Eq. 1)2.Ill) holds, the topological term is 
absent, and the effective action in the partition function 
is simply the action (12.81) for a QNLSM on a spherei^ In 
this case, the quantum-disordered phase at the strong¬ 
coupling fixed point g —^ cx) is stable. In this strongly in¬ 
teracting phase and when © = Z, quantum fluctuations 
restore dynamically and non-perturbatively all the sym¬ 
metries broken by the saddle-point, including any pro¬ 
tecting symmetries. If so, all boundary Dirac fermions 
are gapped out. We then conclude that 


The stability 


®int = ® (2.14) 

when © = Z 2 follows from the fact that one of the ho¬ 
motopy groups with D < d -|- 1 is always 

non-trivial when © = Z 2 (see Sec. 1111 Dl) . 

As an illustration of this method, we give in Table U 
the equivalence classes of topological insulators and su¬ 
perconductors belonging to the ten AZ symmetry classes 
in the presence of interactions that select a short-ranged 
entangled many-body ground state. It becomes apparent 
that the Bott periodicity of the tenfold way, i.e., the peri¬ 
odicity of the (zeroth) homotopy groups of the classifying 
spaces with respect to d, is lost. It also becomes apparent 
that the reduction of the topologically distinct equiva¬ 
lence classes of noninteracting fermions for any given AZ 
symmetry class occurs only in odd dimensions of space. 
Finally, two of the AZ symmetry classes, namely the chi¬ 
ral symmetry classes BDI and CII, have the particularity 
that they may be interpreted either as a superconductor 
or an insulator. Gorrespondingly, the reduction of their 
classification Z —^ for the superconductor interpreta¬ 

tion and Z —> Z„ for the insulator interpretation of these 
symmetry classes obeys 


m = 2n, (class BDI) (2.15a) 

m = n, (class GII) (2.15b) 

when the dimensionality of space is d = 1 mod 4. 
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TABLE 1. The ten Altland-Zirnbauer (AZ) symmetry classes and their topological classification when (i) fermion-fermion 
interactions neither break explicitly their defining symmetries nor spontaneously, (ii) and the many-body ground state is short- 
ranged entangled. Two complex and eight real symmetry classes are characterized by the presence or the absence of time-reversal 
symmetry (T), particle-hole symmetry (C), and chiral symmetry (Ts). Their presence is complemented by the sign multiplying 
the identity in = ±1 or = ±1, and by 1 for Ts. Their absence is indicated by 0. For each symmetry class and for 
any dimension d = 0,1, 2,... of space, the classifying space V^, the space of normalized Dirac masses allowed by symmetry, is 
given in the fifth column. Explicit forms of the classifying spaces C, and and their stable homotopy groups are found in 
Table IXVTl from Appendix]^ The reduction, if any, that arises from the effects of interactions on the topological classification 
of noninteracting fermions for d = 1,...,8 is given in the last eight columns. Each entry with a non-trivial Abelian group 
defines equivalence classes of interacting topological insulators (superconductors) with a short-ranged entangled many-body 
ground state. We color in blue the entry corresponding to a given symmetry class and a given column of odd dimensionality 
d to indicate that this entry is a quotient group of 0 = Z. The reduction Z 0;^^ results from an instability of the 
noninteracting topological classification to fermion-fermion interactions. The four entries corresponding to the symmetry classes 
BDI and CII and the dimensions d = 1 and d — 5 occur in pairs depending on whether these two classes are interpreted as 
describing superconductors (i.e., interacting Majorana fermions) or insulators (i.e., interacting complex fermions), respectively. 


Class 

T 

c 

Ts 


d = 1 

d = 2 

d = 3 

d = 4 

d = 5 

d = 6 

II- 

II 

d = 8 

A 

0 

0 

0 

Co+d 

0 

Z 

0 

z 

0 

Z 

0 

z 

AIII 

0 

0 

1 

Ci+d 

^4 

0 

Zg 

0 

Zig 

0 

Z32 

0 

AI 

+1 

0 

0 

Ro-d 

0 

0 

0 

z 

0 

Z2 

Z2 

z 

BDI 

-fl 

-n 

1 

R-l-d 

Zg, Z4 

0 

0 

0 

^ 16 ! ^8 

0 

Z2 

Z2 

D 

0 

-bl 

0 

R2-d 

Z2 

z 

0 

0 

0 

z 

0 

Z2 

Dill 

-1 

-bl 

1 

Rs-d 

Z2 

Z2 

^16 

0 

0 

0 

^32 

0 

All 

-1 

0 

0 

Ri-d 

0 

Z2 

Z2 

z 

0 

0 

0 

z 

CII 

-1 

-1 

1 

Rs-d 

Z2, Z2 

0 

Z2 

Z2 

Z161 ^16 

0 

0 

0 

C 

0 

-1 

0 

^ 6 —d 

0 

z 

0 

Z2 

Z2 

z 

0 

0 

Cl 

-bl 

-1 

1 

Rj-d 

0 

0 

Z4 

0 

Z2 

Z2 

Z32 

0 


III. REDUCTION OF THE PERIODIC TABLE 
FOR STRONG TI AND TS 


A. The case of one-dimensional space 


In this section, we apply the strategy explained above 
to study the breakdown of the tenfold way in the presence 
of quartic contact interactions in the ascending order of 
the spatial dimension d, i.e., d = 1,2,3, and higher di¬ 
mensions. 

We will use the following conventions. The operation of 
complex conjugation will be denoted by K. Linear maps 
of two-dimensional vector space shall be represented 
by 2 X 2 matrices that we expand in terms of the unit 
matrix Tq and the three Pauli matrices T 2 , and Tg. 
Linear maps of the four-dimensional vector space = 
(8> will be represented by 4 x 4 matrices that we 
expand in terms of the 16 Hermitian matrices 

® V’ = 0,1,2,3, (3.1) 

where cr^, is a second set comprised of the unit matrix 
and the three Pauli matrices. Linear maps of the 2"’- 
dimensional vector space <8) • • • <8> will be 

represented by 2” x 2" matrices that we expand in terms 
of the 4" Hermitian matrices 


A 




Ml 


M2 


where /ii,--- = 0,1, 2, 3. 


(3.2) 


Fidkowski and Kitaev showed in Ref. [13 that, in one 
spatial dimension, any pair of Hamiltonian in the sym¬ 
metry class BDI whose noninteracting topological indices 
differ by eight can be transformed into each other adia- 
batically (i.e., without closing the spectral gap) in the 
presence of a quartic contact interaction that preserves 
TRS. This work was followed up in Refs. [13 and [13 with 
the construction of a topological invariant for interact¬ 
ing fermions from the matrix product representation of 
ground states. This topological invariant establishes that 
the reduction Z —^ Zg is exhaustive. The same approach 
with matrix product states was used to obtain an exhaus¬ 
tive classification of one-dimensional gapped spin systems 
in Ref. P 

Here, we focus on the three chiral symmetry classes 
that support the Z topological classification in the nonin¬ 
teracting limit. We shall reproduce the reduction Z Zg 
and Z ^ Z 2 when the symmetry classes BDI and CII are 
interpreted as chains of Majorana fermions, respectively. 

The one-dimensional chiral symmetry classes can be 
also realized as chains of complex fermions with sub¬ 
lattice symmetry and fermion-number conservation, e.g., 
polyacetylene. For example, polyacetylene-like chains re¬ 
alize the symmetry class AHI when TRS is broken, the 
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symmetry class BDI when both TRS and the SU (2) spin- 
rotation symmetry are present, and the symmetry class 
CII when TRS holds but not the SU{2) spin-rotation 
symmetry if spin-orbit coupling is sizable. We show that 
the reduction of the noninteracting topological classifica¬ 
tion is Z —Z 4 for the symmetry classes AIII and BDI, 
while it is Z Z 2 for the symmetry class CII, provided 
conservation of the fermion number holds. 


1. The symmetry class BDI when d = 1 

Consider the one-dimensional bulk single-particle 
Dirac Hamiltonian (with Dirac matrices of dimension 
r = 2 = r^-J, 

:=-id,T, + m{x)T 2 . (3.3a) 

This single-particle Hamiltonian belongs to the symme¬ 
try class BDI, for 

TU^°\x) = +n^°Hx), (3.3b) 

Cn^°'>{x)C-^ =-H^°\x), (3.3c) 

where 


where 'tisa.vxv unit matrix. Observe that T and C com¬ 
mute with 0 1 and with each other. The domain wall 
(I3.4all must then support v linearly independent bound¬ 
ary zero modes. They are annihilated by the boundary 
Hamiltonian 

= ^bd ® a = 0. (3.6) 

The topological sectors for noninteracting Hamiltonians 
are thus labeled by the integer v taking values in Z in 
the limit v ^ 00 . 

A generic local quartic interaction that respects the 
defining BDI symmetries with the potential to gap out 
these boundary zero modes reduces to a dynamical Dirac 
mass (that depends on imaginary time r in addition 
to space x) that belongs to the symmetry class D, 
upon performing a Hubbard-Stratonovich transforma¬ 
tion. Hence, we must consider the dynamical bulk single¬ 
particle Hamiltonian 

-^Cyn) _|_ rn[x) T 2 ] 01 -I- V(t, x). (3.7a) 

The dynamical Dirac mass V(t, x) is here defined by the 
condition that it anticommutes with H*-°^(x) 0 1 , when 
independent of x, and obeys the transformation laws dic¬ 
tated by the symmetry class D, i.e., it is of the form 


T := Ti K, 


C := To K. (3.3d) 


V{t,x) := 0 7'(t, x), 


7 '(t, x) := iM(T, x), (3.7b) 


The Dirac mass matrix T 2 is here the only one allowed 
for dimension two Dirac matrices under the constraints 
(I3.3bl) and (I3.3cll . As was shown by Jackiw and Rebbi, 
if translation symmetry is broken by the mass term sup¬ 
porting the domain wall 

m(x) = sgn(x), € M, (3.4a) 


where 

M{t,x) = M*{t,x), M{t,x) = —M~^{t,x), (3.7c) 

is a real-valued antisymmetric v x v matrix. Conse¬ 
quently, TRS is only retained for a given V(t, x) if 

M{t,x) = —M{—t,x). (3.7d) 


at X = 0 , then the zero mode 


-i’^3 "^2 / da:' m ( x ') 
e 0 ^ 




(3.4b) 


On the boundary, the operations for reversal of time 
and charge conjugation are now represented by 


7bd •— K, Cdd K. (3.8a) 


where 

HX = sgn(m^)x, (3.4c) 

is the only normalizable state bound to this domain wall. 
This boundary state is a zero mode. It is an eigenstate 
of the single-particle boundary Hamiltonian 

= 0. (3.4d) 


Hence, we must consider the dynamical single-particle 
boundary Hamiltonian 

:=iM(r), (3.8b) 

where M (r) is a real-valued antisymmetric v x v matrix. 
The space of boundary normalized Dirac mass matrices 
obtained by demanding that 7 ' square to the unit v x v 
matrix is topologically equivalent to the space 


Suppose that we consider v = 1 , 2, • • • identical copies 
of the single-particle Hamiltonian (13.31) by defining 

•H(,°)(x) :='H(°)(x )0 1, (3.5a) 

and 

T := T 4 0 1 K, C := Tq 0 1 K, (3.5b) 


K = 0{v)lU{ul2) (3.9) 

for the symmetry class D in zero-dimensional space, pro¬ 
vided the rank u >2 and v is even. The limit ^ ^ 00 of 
these spaces is the classifying space i? 2 - order to gap 
out dynamically the boundary zero modes without break¬ 
ing the defining symmetries of the symmetry class BDI, 
we need to construct a (O-l-I)-dimensional QNLSM for 
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the (boundary) dynamical Dirac masses from the zero¬ 
dimensional symmetry class D without topological ob¬ 
structions. We construct explicitly the spaces for the rel¬ 
evant normalized boundary dynamical Dirac masses of 
dimension ^ = 2" with n = 0,1,2,3 in the followingi^ 
The relevant homotopy groups are given in Table 

Case V = 1: No Dirac mass is allowed on the boundary, 
because the boundary is the end of a one-dimensional 
^2 topological superconductor in the topologically non¬ 
trivial phase of the symmetry class D. 

Case V = 2: We use the representation 1 = ctq. There 
is one dynamical normalized Dirac mass on the boundary 
that is proportional to the matrix (72 ■ ^ domain wall 
in imaginary time such as m 2 oo sign(T) (72 prevents the 
dynamical generation of a spectral gap on the boundary. 

Case V = A: We use the representation 1 = CTq (g) Pq. 
A (maximum) set of pairwise anticommuting boundary 
dynamical Dirac mass matrices follows from the set 

{(T2® P2}- ( 3 - 10 ) 

This set spans the space of normalized boundary dynam¬ 
ical Dirac masses that is homeomorphic to . Even 
though 7rg^2(‘5'^) = 0; it is possible to add a topolog¬ 
ical term that is nonlocal, yet only modifies the equa¬ 
tions of motion of the (O-l-l)-dimensional QNLSM on the 
boundary by local terms as a consequence of the fact that 
7 rQ+i+i(5'^) = Z. Such a term is a (O-l-l)-dimensional ex¬ 
ample of a Wess-Zumino (WZ) term. In the presence 
of this WZ term, the boundary theory remains gapless. 
It is nothing but a bosonic representation of the gapless 
S' = 1 /2 degrees of freedom at the end of a quantum spin- 
1 antiferromagnetic spin chain in the Haldane phasei^^ 
Case V = 8.' We use the representation 1 = Uq^Pq^Aq. 
One set of pairwise anticommuting boundary dynamical 
Dirac mass matrices follows from the set 

{(72 ® Pq O Ao,( 73 (g)p2 ® Ao,( 73 (8IP3 O A2,(73 (g) Pg (g> A 2 } 

( 3 . 11 ) 


TABLE II. Reduction from Z to for the topologically 
equivalent classes of the one-dimensional SPT phases in the 
symmetry class BDI that arises from interactions. We de¬ 
note by W the space of v x v normalized Dirac mass matrices 
in zero-dimensional Hamiltonians belonging to the symmetry 
class D. The limit v ^ 00 of these spaces is the classifying 
space R 2 - The second column shows the stable D-th homo¬ 
topy groups of the classifying space R 2 . The third column 
gives the number u of copies of boundary (Dirac) fermions 
for which a topological obstruction is permissible. The fourth 
column gives the type of topological obstruction that prevents 
the gapping of the boundary (Dirac) fermions. 


D 

■^£1(7^2) 

V 

Topological obstruction 

0 

Z2 

2 

Domain wall 

1 

0 



2 

z 

4 

WZ term 

3 

0 



4 

0 



5 

0 



6 

z 

8 

None 

7 

Z2 




2. The symmetry class CII when d = 1 

Consider the one-dimensional bulk single-particle 
Dirac Hamiltonian (with Dirac matrices of dimension 
r = A = r^i„), 

:= -id^ A 30 -b m{x) X^q. (3.12a) 

This single-particle Hamiltonian belongs to the symme¬ 
try class CII, for 

T-H^°\x) r-^ = +n^°Hx), (3.12b) 

cn^°'>{x)c-^ = ( 3 . 12 c) 

where 

r := 1 X 32 K, C := iAg 2 K. (3.12d) 


This set spans a manifold homeomorphic to (we may 
find a set of pairwise anticommuting masses spanning 
5®). No topological term is admissible over this target 
manifold that delivers local equations of motion. The 
QNLSM over this target space endows dynamically the 
boundary Hamiltonian with a spectral gap. 

We conclude that the effects of interactions on the one¬ 
dimensional SPT phases in the symmetry class BDI are 
to reduce the topological classification Z in the nonin¬ 
teracting limit down to Zg under the assumption that a 
Hamiltonian from the symmetry class BDI is interpreted 
as a mean-field description of a superconductor. The 
logic used to reach this conclusion is summarized by Ta¬ 
ble H once the line corresponding to ix = 2 has been 
identified. It is given by the smallest D that accommo¬ 
dates a non-trivial entry for the corresponding homotopy 
group. The line for zz = 4 is then identified with the next 
smallest D with 7r^(i?2) ^ 0; so on. 


The Dirac mass matrix A 2 g is here the only one allowed 
for dimension four Dirac matrices in the symmetry class 
CII. If translation symmetry is broken by the Dirac mass 
term supporting the domain wall 


"^(2;) ='^ooSgn(a;), 

at a; = 0, then the zero mode 


-i-^30 Wo / da;'m(a;') 

e » y = 


€ M, 


x\ 


X, 


where 


(3.13a) 


(3.13b) 


^10 X = sgn (m,^) X, (3.13c) 

is the only normalizable state bound to this domain wall. 
This boundary state is a zero mode. It is an eigenstate 
of the single-particle boundary Hamiltonian 


(3.13d) 







Suppose that we consider v = 1, 2, • • ■ identical copies 
of the single-particle Hamiltonian (13.31) by defining 

(3.14a) 

and 

r := iXi2 0 1 K, C := iXo2 0 1 K, (3.14b) 

where H is a v x v unit matrix. Observe that 'T and C 
commute with X^q 0 1 and with each other. The do¬ 
main wall (I3.13a|) must then support v linearly indepen¬ 
dent boundary zero modes. They are annihilated by the 
boundary Hamiltonian 

® 1 = 0 - ( 3 - 15 ) 

The topological sectors for noninteracting Hamiltonians 
are thus labeled by the integer v taking values in Z in 
the limit v ^ oo. 

A generic local quartic interaction that respects the 
defining CH symmetries with the potential to gap out the 
boundary zero modes reduces to a dynamical Dirac mass 
(that depends on imaginary time r in addition to space x) 
that belongs to the symmetry class C, upon performing a 
Hubbard-Stratonovich transformation. Hence, we must 
consider the dynamical bulk single-particle Hamiltonian 

(t, x) := [-id^ Ago -h m{x) A 20 ] 0 1 -b V(r, x), 

(3.16a) 

where the dynamical Dirac mass V(t, x) is defined by the 
condition that it anticommutes with TL^^^x) 0 1, when 
independent of x, and obeys the transformation laws dic¬ 
tated by the symmetry class C, i.e., it must obey 

C V(t, x) = — V(r, x). (3.16b) 

On the boundary, the operations for reversal of time 
and charge conjugation are now represented by 

7^0 := 0-2 0 1 K, Cod := 0-2 0 11 K. (3.17a) 

Hence, we must consider the dynamical single-particle 
boundary Hamiltonian 

:= 7'(t), (3.17b) 

where 

Cbd7'('r)Cod = -I'ir). (3.17c) 

The space of normalized Dirac mass matrices obtained by 
demanding that 7 ^(r) squares to the unit 2iy x 2u matrix 
for all imaginary times is the space 

K = Sp{y)/U{u) (3.18) 

for the symmetry class C in zero-dimensional space. The 
limit > 00 of these spaces is the classifying space i?g. 
In order to gap out dynamically the boundary zero modes 


TABLE III. Reduction from Z to Z 2 for the topologically 
equivalent classes of the one-dimensional SPT phases in the 
symmetry class CII that arises from interactions. We denote 
by the space oi v x v normalized Dirac mass matrices 
in zero-dimensional Hamiltonians belonging to the symmetry 
class C. The limit zz —>■ 00 of these spaces is the classifying 
spaces Rq. The second column shows the stable D-th homo- 
topy groups of the classifying space Ag- The third column 
gives the number v of copies of boundary (Dirac) fermions 
for which a topological obstruction is permissible. The fourth 
column gives the type of topological obstruction that prevents 
the gapping of the boundary (Dirac) fermions. 


D 

'^oiRe) 

V 

Topological obstruction 

0 

0 



1 

0 



2 

z 

1 

WZ term 

3 

Z 2 

2 

None 

4 

Z 2 



5 

0 



6 

z 



7 

0 




without breaking the defining symmetries of the symme¬ 
try class CH, we need to construct a (O-l-l)-dimensional 
QNLSM for the (boundary) dynamical Dirac masses from 
the zero-dimensional symmetry class C without topo¬ 
logical obstructions. We construct explicitly the spaces 
for the relevant normalized boundary dynamical Dirac 
masses of dimension zz = 2” with n = 0,1 in the follow¬ 
ing. The relevant homotopy groups are given in Table iHll 

Case V = \: The three Dirac mass matrices cr^, (T 2 , and 
CTg, are allowed on the boundary. They all anticommute 
pairwise. A WZ term is permissible as (5^) = Z. 

In the presence of this WZ term, the boundary theory 
remains gapless. 

Case V = 2: The minimum number of anticommuting 
mass matrices is larger than three. Hence, the zeroth, 
first, and second homotopy groups over the boundary 
normalized dynamical Dirac masses all vanish. No topo¬ 
logical term is possible. The (O-l-l)-dimensional QNLSM 
over this target space endows dynamically the boundary 
Hamiltonian with a spectral gap. 

We conclude that the effects of interactions on the one¬ 
dimensional SPT phases in the symmetry class CH are 
to reduce the topological classification Z in the nonin¬ 
teracting limit down to Z 2 under the assumption that a 
Hamiltonian from the symmetry class CH is interpreted 
as a mean-field description of a superconductor. The 
logic used to reach this conclusion is summarized by Ta¬ 
ble [m] once the line corresponding to zz = 1 has been 
identihed. It is given by the smallest D that accommo¬ 
dates a non-trivial entry for the corresponding homotopy 
group. 
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3. The chiral symmetry classes as one-dimensional 
insulators 


from a local quartic interactions after a Hubbard- 
Stratonovich transformation takes the form (I3.7al) with 


So far we have interpreted the symmetry classes BDI 
and CII as examples of topological superconductors by 
focusing on the fact that their second-quantized Hamil¬ 
tonian respects a unitary charge-conjugation symmetry, a 
PHS (see Appendix]^. As the symmetry classes BDI and 
CII also preserve TRS, the composition of reversal of time 
with charge conjugation delivers a non-unitary symme¬ 
try of their second-quantized Hamiltonian, namely a CHS 
(see Appendix lAT) . The third chiral symmetry class AHI 
is defined by demanding that it preserves CHS, no more 
and no less. Hence, any representative gapped Hamilto¬ 
nian from the chiral symmetry class AHI can always be 
interpreted as a topological insulator with fermion num¬ 
ber conservation. 

The CHS can be implemented at the single-particle 
level by a unitary sublattice spectral symmetry for (com¬ 
plex) electrons hopping between two sublattices. From 
this point of view, the three chiral symmetry classes AHI, 
BDI, and CII, when interpreted as metals or as insula¬ 
tors, can be treated on equal footing. The fermion num¬ 
ber is conserved in a metal or in an insulator, unlike in 
the mean-field treatment of a superconductor. This is the 
case when the symmetry class BDI is interpreted as an 
effective theory for polyacetylene, in which case the Dirac 
gap is induced by coupling the electrons to phonons, i.e., 
it realizes a Peierls or bond-density wave instabilityi^ In 
this interpretation of the chiral classes, it is necessary 
to introduce an additional particle-hole grading in order 
to include through dynamical Dirac masses the effects 
of superconducting fluctuations induced by any quartic 
interaction. Failure to do so can produce a distinct re¬ 
duction pattern of the noninteracting topological equiv¬ 
alence classes arising from interactions, for it can matter 
whether the boundary dynamical masses belong to the 
classifying space associated with the symmetry classes D 
or to the symmetry class A. 

a. Symmetry class AIII Consider the one¬ 
dimensional bulk single-particle Dirac Hamiltonian 
in the symmetry class AHI 

:= -ia.Tg -\-m{x)T 2 . (3.19a) 

It anticommutes with the unitary operator 

r 5 :=Ti. (3.19b) 

It supports the zero mode (13.41) at the boundary where 
it identically vanishes, 


V{t,x) := Ti j'{t,x) (3.20) 


and 7 '(t, x) a, v x v Hermitian matrix. 

On the boundary, we must consider the dynamical 
single-particle boundary Hamiltonian 

^m“^M=7'(t). (3.21) 

The u X u Hermitian matrix 7 '(t) belongs to the zero¬ 
dimensional symmetry class A. Consequently, it is as¬ 
signed the classifying space Cg in the limit u ^ oo with 
the zeroth-homotopy group t^q{Cq) = Z. When v = 

7' is a real number, and the domain wall 7^(t) oc sgn(r) 
binds a zero mode at r = 0 [i.e., a normalizable zero 
mode of the operator When v = 2, 

we can write 7 ^(t) as a linear combination of the ma¬ 
trices (7^ with the real-valued functions rn^ir) as coef¬ 
ficients for fx = 0,1, 2,3, respectively. Any one of the 
three Pauli matrices (o’j^, cr 2 , Ug) anticommutes with the 
other two Pauli matrices. Hence, the space of normal¬ 
ized boundary dynamical Dirac masses that anticommute 
pairwise is homeomorphic to with the homotopy group 
7 ro_|_i_|_i( 5 '^) = Z when v = 2. A (O-l-l)-dimensional 
QNLSM for the (boundary) dynamical Dirac masses is 
augmented by a WZ term. 

We conclude that the effects of interactions on the one¬ 
dimensional SPT phases in the symmetry class AHI are 
to reduce the topological classification Z in the noninter¬ 
acting limit down to Z 4 under the assumption that only 
fermion-number-conserving dynamical Dirac masses are 
included in the single-particle Hamiltonian. The logic 
used to reach this conclusion is summarized by Table IIVI 
once the line corresponding Xo v = \ has been identi¬ 
fied. It is given by the smallest D that accommodates a 
non-trivial entry for the corresponding homotopy group. 

To include the effects of superconducting fluctua¬ 
tions in the single-particle Hamiltonian after a Hubbard- 
Stratonovich transformation, we need to consider the di¬ 
rect sum 


H^°j^(x) := H(°)(x)(8)]1 © -H(°)*(x)®]1 


= n^^\x) ® t ® Pq. 


(3.22a) 


This Bogoliubov-de-Gennes (BdG) single-particle Hamil¬ 
tonian anticommutes with the operation for charge con¬ 
jugation 


C := Tg © 1 © P 4 K, (3.22b) 


H[,°i(x)=0. (3.19c) 

Upon tensoring Hamiltonians (I3.19al) and (I3.19cl) to¬ 
gether with the Dirac Fg matrix by the v xv unit matrix 
1, there follows ;/ = 1, 2, 3, • ■ • bonndary zero modes. 

The dynamical single-particle Hamiltonian that en¬ 
codes those non-superconducting fluctuations arising 


in addition to anticommuting with © 1 © pg. Hence, 
it belongs to the symmetry class BDI. The dynamical 
single-particle Hamiltonian that accounts for supercon¬ 
ducting fluctuations (dynamical Dirac masses from the 
symmetry class D) takes the form 

V(r, x) := © 7 '(t, x). 


(3.23a) 
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TABLE IV. Reduction from Z to Z 4 for the topologically 
equivalent classes of the one-dimensional SPT phases in sym¬ 
metry classes AIII or BDI that arises from the fermion- 
number-conserving interacting channels. We denote by 
the space oi v x u normalized Dirac mass matrices in zero¬ 
dimensional Hamiltonians belonging to the symmetry class A. 
The limit iz —>■ 00 of these spaces is the classifying spaces Cq. 
The second column shows the stable D-th homotopy groups of 
the classifying space Cq. The third column gives the number 
v of copies of boundary (Dirac) fermions for which a topolog¬ 
ical obstruction is permissible. The fourth column gives the 
type of topological obstruction that prevents the gapping of 
the boundary (Dirac) fermions. 


D 

'^d(C'o) 

V 

Topological obstruction 

0 

Z 

1 

Domain wall 

1 

0 



2 

z 

2 

WZ term 

3 

0 



4 

z 

4 

None 

5 

0 



6 

z 



7 

0 




where the 2i> x 2v Hermitian matrix must obey 

7'(t, x) = -(tq (g) ]1(8)/9 i)7'*(t, x) (tq (g) lligjpi). (3.23b) 

The stability analysis of the boundary zero modes is sim¬ 
ilar to the one performed below Eq. (I3.9|l except that one 
must replace v in Eq. (13.91) by = 2v and that we 
have a different representation of the PHS. When jz = 1, 

i'{t) = M{t) M(r)e]R, (3.24) 

supports a domain wall in imaginary time on the bound¬ 
ary. When ^ = 2, we use the representation 1 = ctq 
and introduce the notation p^, with p, p' = 

0,1, 2, 3. Now, 

7 '(r)= ^ M^^,(r)eK, (3.25) 

where the sum on the right-hand side is to be per¬ 
formed over the six matrices X 2 i,X 22 ,Xq^,Xi^,X 2 q,X^^. 
This set of six matrices decomposes into two triplets 
of pairwise anticommuting matrices. The first triplet 
is given by X 2 i,X 22 ,Xq^. The second triplet is given 
by Xi^,X 2 Q,X^^. Each triplet defines a two-sphere . 
Hence, each triplet has the potential to accommodate 
a WZ term. However, we must make sure that the in¬ 
tegrity of any one S"^ entering the decomposition 5^ U 5^ 
of the normalized dynamical masses on the boundary is 
compatible with maintaining the global 17(1) symmetry 
associated with the conservation of the fermion number. 

To address the fate of the fermion-number conserva¬ 
tion, observe that Hamiltonian (I3.19al) is invariant under 
the global 17(1) transformation 

17(a) := (3.26) 


with 0 < a < 27r independent of x. This symmetry is 
to be preserved when treating superconducting fluctua¬ 
tions. In the BdG representation (l3.22aL this symmetry 
becomes the symmetry under the global 17(1) transfor¬ 
mation 

<dG ^ Z^BdG(«) ^BdG ^B'dG(«)> (3-27a) 

where 

WBdG(a) := (3.27b) 

When zz = 1, the boundary dynamical mass p.24l) 
is invariant under multiplication from the left with 
exp(-(-ia P 3 ) and multiplication from the right with 
exp(—iapg). When zz = 2, if we normalize the bound¬ 
ary dynamical mass (13.251) by demanding that 

1 = M|i + M |2 + M 23 , (3.28a) 

1 = M23 + a4 + a 4, (3.28b) 

we may then identify these boundary dynamical masses 
as the union of two two-spheres . The global 17(1) 
transformation defined by multiplication from the left 
with exp(-|-ia Xq^) and multiplication from the right with 
exp(—iaXpg) leaves the two-sphere (I3.28al) invariant as 
a set, for it is represented by a rotation about the north 
pole Xgg that rotates the equator spanned by X 21 and 
X 22 with the angle 2a. The same transformation leaves 
the two-sphere (I3.28b|) invariant point-wise. Hence, each 
in U is compatible with the conservation of the 
global fermion number. Because 7 rQ_|_ 3 _|_ 3 (S'^) = Z, a WZ 
topological term in the QNLSM for the boundary is per¬ 
missible. 

We may then safely conclude that the effects of inter¬ 
actions on the one-dimensional SPT phases in the sym¬ 
metry class AIII are also to reduce the topological clas¬ 
sification Z in the noninteracting limit down to Z 4 under 
the assumption that superconducting fluctuation chan¬ 
nels are included in the stability analysis. The Z 4 classi¬ 
fication for one-dimensional SPT phases in the symmetry 
class AIII agrees with the one derived using group coho¬ 
mology in Ref. [with AIII interpreted in Ref. | 6 l| as 
a time-reversal-symmetric superconductor with the full 
spin- 1/2 rotation symmetry broken down to a C/(l) sub¬ 
group]. 

b. Symmetry class BDI Dirac Hamiltonians in the 
symmetry class BDI are obtained from those in the sym¬ 
metry class AIII by imposing the constraint of TRS, Eqs. 
(I3.3bl) and (I3.3dl) (note that = 2 for both AIII and 
BDI classes in d = 1). Since the TRS is not relevant 
for dynamical Dirac masses in the single-particle Dirac 
Hamiltonian after a Hubbard-Stratonovich transforma¬ 
tion, the stability analysis of gapless boundary states 
in the symmetry class BDI in d = 1 follows from that 
of the symmetry class AIII. Consequently, the effects 
of interactions in the symmetry class BDI, when inter¬ 
preted as realizing complex fermions as opposed to Ma- 
jorana fermions, is to reduce the topological classification 
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Z in the noninteracting limit down to Z 4 under the as¬ 
sumption that only fermion-number preserving dynami¬ 
cal Dirac masses taken from the symmetry class A are 
included in the stability analysis. Furthermore, if dy¬ 
namical superconducting fluctuations are allowed by in¬ 
troducing an additional particle-hole grading and dynam¬ 
ical Dirac masses from the symmetry class D, then the 
same reduction pattern Z —>■ Z 4 follows. The topological 
classification Zg of the symmetry class BDI when inter¬ 
preted as describing Majorana fermions is thus finer than 
the classification Z 4 of the symmetry class BDI when in¬ 
terpreted as describing complex fermions. 

c. Symmetry class CII We interpret the single¬ 
particle Hamiltonian ()3.12a|) as describing an insulator, 
not a superconductor. This is to say that the defining 
symmetries are TRS (I3.12bl) and the CHS 

Tg r-i = Tg := Alo. (3.29) 

If we are after the dynamical effects of interactions that 
preserve the (complex) fermion number, we may use Eq. 
(13.161) with the only caveat that the dynamical mass ma¬ 
trix is now required to belong to the symmetry class A in¬ 
stead of the symmetry class C. The boundary dynamical 
Hamiltonian is then the same as for the symmetry classes 
AHI and BDI, i.e., Eq. (13.211) . except for its rank being 
twice as large as compared to the symmetry classes AHI 
and BDI, for an additional grading (that of the spin-1/2 
degrees of freedom) has been accounted for. This larger 
rank implies that the WZ term is already permissible at 
V = 1, i.e., the reduction pattern for the noninteracting 
topological classihcation is Z ^ Z 2 . It remains to verify 
that the same reduction pattern is also obtained if su¬ 
perconducting fluctuations are included, as was the case 
for the symmetry classes AHI and BDI. To this end, we 
must tensor (|3.12aD with pg, in which case the charge 
conjugation symmetry is realized by Tq ® ctq 0 Pi K. We 
can borrow the stability analysis with respect to super¬ 
conducting interacting channels from the symmetry class 
D that we performed for the symmetry classes AHI and 
BDI, again with the caveat that the rank of the boundary 
dynamical Hamiltonian is twice as large as it was. This 
larger rank implies again that the WZ term is already 
permissible for v = 1 , i.e., the reduction pattern for the 
noninteracting topological classification is again Z Z 2 . 
Thus we obtain the same topological classification Z 2 of 
the symmetry class CII in d = 1 both when interpreted 
as describing Majorana fermions (superconductors) and 
when interpreted as describing complex fermions (insu¬ 
lators). 


B. The case of two-dimensional space 


to average the Kubo Hall conductivity over all twisted 
boundary conditions of the many-body ground state as 
a signature of both the IQHE and EQHE. A mathemati¬ 
cally ri gor ous proof of this intuition can be found in Refs. 

andl63l. This intuition readily extends to the symme¬ 
try class D and C as they realize quantized thermal Hall 
effects. The robustness of chiral edge modes in the sym¬ 
metry classes D, C, and A to quartic contact interactions 
will be derived using the method of Sec. HIl 

Let X = (xi,X 2 ) denote a point in two-dimensional 
space. The single-particle Dirac Hamiltonian with the 
smallest rank = 2 that admits a Dirac mass can be 
chosen to be represented by 

^A^(a:) := + ^i(a:)] ^3 -k [-i^a + A^ix)] 

+ Aq{x) + m{x) (72- (3.30) 


It belongs to the symmetry class A for arbitrarily chosen 
vector potentials A(a;), scalar potential Aq{x), and mass 
m{x). When the gauge fields are vanishing. 


:= -i^iCTa 


\d 2 Oi + m{x) (T 2 


n^°\x) 


(3.31) 


belongs to the symmetry class D. Finally, the single¬ 
particle Hamiltonian with the smallest rank = 4 
that belongs to the symmetry class C can be chosen to 
be represented by 



3 

-i^iAao-f ^Ay(a;)A3^- 

i=i 

3 

-i92^io + X!^2j(a:)Ay 

i=i 

3 

+ X! ^Oj + Mx) X20 


-A, 


02 


n^°\x) 


A, 


02 - 


(3.32) 


In two spatial dimensions the symmetry classes A, 
D, and C realize noninteracting topological insula¬ 
tors and superconductors with the Grassmannian man¬ 
ifolds Cg = limAr_),oo u()Lo 17(A)/[C/(n) x U{N - n)], 
i?o = limAr^ooU))LgO(A)/[ 0 (n)x)(A - n)], and R 4 = 
limjv->.oo U))Lq «S'p(A)/[S'p(n)x)(A — n)] as classifying 
spaces, respectively. They share the same zeroth-order 
homotopy group Z. This group also serves as defining 
the topological attributes of noninteracting topological 
insulators and superconductors in the symmetry classes 
A, D, and C. 


The notion that the chiral edge modes in the IQHE are 
immune to local interactions is rather intuitive. Neither 
backscattering nor umklapp scattering is permissible. An 
operational and quantitative validation for this intuition 
goes back to Niu and Thouless in Ref. |^, who proposed 


1. The symmetry class D when d = 2 

Let 1 denote a u x v unit matrix with v = 1,2, • • •. 
Consider the two-dimensional bulk single-particle Dirac 
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Hamiltonian 

T-L^^^x) := —i9;^ Cg®! — ii 92 cr^ ( 8 ) 1+ m(a:) 0-2 (3.33) 

of rank 2 v. There is no Hermitian (2z/) x { 2 v) matrix 
that anticommutes with If so, the set {/3} in 

Eq. (12.31) is empty. In other words, no dynamical mass 
is available to induce a dynamical instability of the v 
boundary zero modes. 


2. The symmetry class C when d = 2 


The same reasoning applies to the bulk single-particle 
Hamiltonian 


(x) := [-i^i X30 - ia2 + m{x) X^q] 1, 

(3.34a) 

of rank 4zz that realizes a topological superconductor in 
the symmetry class C, 


^ ^ r^( 0 ) ( 3 .)! ^Xo 2 ® a) • (3.34b) 


No dynamical mass is available to induce a dynamical 
instability of the ly boundary zero modes. 


3. The symmetry class A when d — 2 

If the single-particle Dirac Hamiltonian (13.331) is inter¬ 
preted as describing an insulator with fermion-number 
conservation, then no dynamical mass that preserves the 
fermion number and anticommutes with (T 2 ® 1 is per¬ 
missible. The same remains true if we account for super¬ 
conducting fluctuations, for the BdG extension of (13.331) 
that is given by Eq. (I3.34al) . whereby charge conjugation 
is defined by [and not by Eq. (I3.34b|) ] 

C:=Xoi(8)1K, (3.35) 

fails to anticommute with any (4j^) x (4:/) Hermitian ma¬ 
trix allowed by the PHS generated by the operation of 
charge conjugation (13.351) . 


C. The case of three-dimensional space 

The reduction Z —>■ Z^g for the three-dimensional in¬ 
teracting topological superconductors belonging to the 
symmetry class Dili has been understood in the fol¬ 
lowing ways after a conjecture by Kitaev from Ref. [isl . 
One approach is to enumerate the distinct topologi¬ 
cal orders at the two-dimensional surface of the three- 
dimensional bulk that cannot be realized with any bulk 
two-dimensional Hamiltonianj^^— In this approach, the 
breakdown of Z takes place when vortices (point-like de¬ 
fects of a symmetry-broken phase) at the surface prolif¬ 
erate (deconhne) so as to stabilize a gapped and fully 


symmetric surface phase. Another approach advocated 
by You and Xu consists in relating fermionic short-ranged 
entangled ground states to bosonic short-ranged entan¬ 
gled ground states 1 ^ They also applied their approach 
to systems with inversion symmetry. The reductions 
Z Z 4 and Z Zg for the symmetry classes Cl and 
AHI was obtained in Ref. [U Finally, Kapustin has pro¬ 
posed to classify symmetry protected topological phases 
for interacting bosons or fermions by considering low- 
energy effective actions that are invariant under cobor- 
dism (a certain type of equivalence relation between man¬ 
ifolds) i^“— 

We apply the method of Sec.|TTlto the symmetry classes 
Dili, Cl, and AHI in the presence of quartic contact in¬ 
teractions. We recover the reductions Z —>• Zgg, Z —>• Z 4 , 
and Z ^ Zg for the symmetry class Dili, Cl, and AHI, 
respectively. We also verify that the topological classifi¬ 
cation Z 2 of the symmetry class AH is stable to quartic 
contact interactions. 

We shall denote with x = {x,y,z) = (xi, X 2 , x^) a 
point in three-dimensional space. 


1. The symmetry class Dill when d = 3 

Let X^^, = with = 0,1,2,3, Consider the 

three-dimensional bulk single-particle Dirac Hamiltonian 
(with Dirac matrices of dimension r = 4 = r’jjjjjj), 

TL^^\x) := —i 9 ]^ X^^ — 182 Xq 2 — id^ A44 -|- •m(x) X^^. 

(3.36a) 

This single-particle Hamiltonian belongs to the three- 
dimensional symmetry class Dili, for 

Tn^°Hx) (3.36b) 

C'H^°\x)C-^ =-VS°\x), (3.36c) 

where 

r := iX 2 o K, C := Agg K. (3.36d) 

The multiplicative factor i in the definition of 'T is needed 
for T to commute with C. 

The Dirac mass matrix Agg is here the only one allowed 
for dimension four Dirac matrices under the constraints 
(I3.36bl) and (I3.36cl) . Consequently, the domain wall 

m{x) = m{y) := sgn(?/), e R, (3.37a) 


at y = 

0 , binds the zero mode 



y 

-iW2-^03 /dy'™(y') . 

e “ X = e- ' 

(3.37b) 

where 

^oiX = -sgn(m^)x 

(3.37c) 


with X independent of x and 2 . The kinetics of the gap¬ 
less boundary states is governed by the Dirac Hamilto¬ 
nian 


(3.38) 
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where we have chosen < 0. 

On the boundary, the operations for reversal of time 
and charge conjugation are now represented by 

1 / := i'r 2 ® 1 K, Cbd 1 . := To (g) 1 K, (3.39a) 

where 1 is the v x v unit matrix. We seek the single¬ 
particle Hamiltonian on the boundary that encodes the 
fluctuations arising from the Hubbard-Stratonovich de¬ 
coupling of quartic interactions through a generic dy¬ 
namical mass that respects the PHS on the boundary. 
It is given by 

X, z) := - id^T3 (g 1 - iS^Ti g 1 

+ T2^M{t,x,z) (3.39b) 

with the v X V real-valued and symmetric matrix 

M(r, X, z) = M*{t, X, z) = M~^{t, x, z). (3.39c) 

The space of normalized Dirac mass matrices of the 
form (I3.39 cI) is topologically equivalent to the space 

1/ 

K := U 0{v)/ [0{k) X 0{v - k)] (3.40) 

for the symmetry class D in two-dimensional space. The 
limit > oo of these spaces is the classifying space Rq. 
In order to gap out dynamically the boundary zero modes 
without breaking the defining symmetries of the symme¬ 
try class Dili, we need to construct a (2-|-l)-dimensional 
QNLSM for the (boundary) dynamical Dirac masses from 
the two-dimensional symmetry class D without topolog¬ 
ical obstructions. We construct explicitly the spaces for 
the relevant normalized boundary dynamical Dirac mass 
matrices of dimension i/ = 2"^ with n = 0,1, 2, 3,4 in the 
following. The relevant homotopy groups are given in 
Table |Vl 

Case ly = 1: There is one boundary dynamical Dirac 
mass matrix ^'{t,x,z) on the boundary that is propor¬ 
tional to r 2 . A domain wall in imaginary time such as 
oo sign(T) T 2 prevents the dynamical generation of a 
spectral gap on the boundary. 

Case V = 2: We use the representation 1 = ctq. The 
2x2 real-valued and symmetric matrix M (r, x, z) is a lin¬ 
ear combination with real-valued coefficients of the pair 
of anticommuting matrices and If M{t,x,z) is 
normalized by demanding that it squares to Cg, then the 
set spanned by M{t,x,z) is homeomorphic to the one- 
sphere S^. As 7r^(S'^) = Z, it follows that M{t,x,z) 
supports vortex lines that bind zero modes in (2-1-1)- 
dimensional space and time and thus prevent the gapping 
of the boundary statesi^ 

Case V = 4; We use the representation 1 = (Tq ® u'q. 

The 4x4 real-valued and symmetric matrix M{t,x,z) 

is a linear combination with real coefficients of A , = 

cr^, 

g a'^, with /i,^' = 0,1,2,3 such that either none or 
two of ^ and /i' equal the number 2. Of these the three 


TABLE V. Reduction from Z to Zig for the topologically 
equivalent classes of the three-dimensional SPT phases in the 
symmetry class Dili that arises from interactions. We de¬ 
note by W the space oi u x u normalized Dirac mass matrices 
in boundary (d = 2) Dirac Hamiltonians belonging to the 
symmetry class D. The limit oo of these spaces is the 
classifying space Rq. The second column shows the stable D- 
th homotopy groups of the classifying space Rq. The third 
column gives the number i/ of copies of boundary (Dirac) 
fermions for which a topological obstruction is permissible. 
The fourth column gives the type of topological obstruction 
that prevents the gapping of the boundary (Dirac) fermions. 


D 

td(Ao) 

V 

Topological obstruction 

0 

Z 

1 

Domain wall 

1 

Z2 

2 

Vortex line 

2 

Z2 

4 

Monopole 

3 

0 



4 

z 

8 

WZ term 

5 

0 



6 

0 



7 

0 



8 

z 

16 

None 


matrices A 33 , and Xqi anticommute pairwise. If 
M(t,x,z) is a linear combinations with real-valued co¬ 
efficients of these three matrices and if M is normal¬ 
ized by demanding that it squares to Agg, then the set 
spanned by M (t, x, z) is homeomorphic to the two-sphere 
S^. As ^ 2 ( 5 ^) = Z, M(t, X, z) supports point-like defects 
of the monopole type that bind zero modes in ( 2 - 1 - 1 )- 
dimensional space and time and thus prevent the gapping 
of the boundary states.— 

Case V = 8 .- We use the representation 1 = o-p g 
CTg g (Tg. The 8 x 8 real-valued and symmetric matrix 
M{t,x,z) is a linear combination with real-valued coef¬ 
ficients of the matrices X g a' , g cr"„ where 

flfl H H fJs 

either none or two of p,', = 0,1, 2, 3 equal the num¬ 

ber 2. Of these, one finds the five pairwise anticom¬ 
muting matrices A 333 , Aj^gg, Agj^g, Aggj^, and X2i2- If 
M{t,x,z) is a linear combination with real-valued coef¬ 
ficients of these five matrices and if M is normalized by 
demanding that it squares to Aggg, then the set spanned 
by M(t, X, z) is homeomorphic to the four-sphere 5^. As 
7 r 4 ( 5 '^) = Z, it is possible to add a topological term to 
the QNLSM on the boundary that is of the WZ type. 
This term is conjectured to prevent the gapping of the 
boundary states. 

Case V = 16.- We use the representation 1 = ctq gagg 
CTg gcTg'. The 16 X 16 real-valued and symmetric matrix 
M(t, X, z) is a linear combination with real-valued coeffi¬ 
cients of the matrices A„ = cr,, g a', g a''„ g 

l-LI-L fJ.fl fl fl fl fl 

where none, two, or four of p,', /i", /r"' = 0 , 1 , 2 ,3 equal 
the number 2. Of these, one finds the nine pairwise anti¬ 
commuting matrices A 22221 ^01221 ^03221 ^20121 ^2032 j 
"^ 1202 ; "^32021 "^oooii ^0003' 1^ M{t,x,z) is a lin¬ 
ear combination with real-valued coefficients of these 
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nine matrices and if M(t, x, z) is normalized by demand¬ 
ing that it squares to ^qqqq, then the set spanned by 
M(t,x,z) is homeomorphic to the eight-sphere 5'®. It is 
then impossible to add a topological term to the QNLSM 
on the boundary. The boundary zero modes can be 
gapped out. 

We conclude that the effects of interactions on the 
three-dimensional SPT phases in the symmetry class Dill 
are to reduce the topological classihcation Z in the non¬ 
interacting limit down to Zj^g. The logic used to reach 
this conclusion is summarized by Table |V] once the line 
corresponding to = 1 has been identified. It is given 
by the smallest D that accommodates a non-trivial en¬ 
try for the corresponding homotopy group. The line for 
ly = 2 is then identified with the next smallest D with 
7r^(i?2) 7^ 0, and so on. 


2. The symmetry class Cl when d = 3 


Let ® with p,v,\ = 0,1,2,3. 

Consider the three-dimensional bulk single-particle Dirac 
Hamiltonian (with Dirac matrices of dimension r = 8 = 

T ■ ) 
min/’ 

(x) := —idiXgiQ — 182 Xq2q — idg X^q + m{x) ^oao- 

(3.41a) 

This single-particle Hamiltonian belongs to the three- 
dimensional symmetry class Cl, for 

rn^^Hx) = +n^°'>{x), (3.41b) 

cn^°'>{x)c-^ =-n^°Hx), (3.41c) 

where 

r := X202 K, C := iXoi2 K. (3.41d) 

The multiplicative factor i in the definition of C is needed 
for T to commute with C. 

The single-particle Hamiltonian p.41al) is the direct 
product of the single-particle Hamiltonian (I3.36al) with 
the unit 2x2 matrix Cg. If we interpret the degrees of 
freedom encoded by Cg and the Pauli matrices cr as car¬ 
rying spin-1/2 degrees of freedom, we may then interpret 
Eqs. (13.411) as defining a spin-singlet superconductor that 
preserves TRS. 

The Dirac mass matrix Wggg is here the only one al¬ 
lowed for dimension eight Dirac matrices under the con¬ 
straints (I3.41bl) and (I3.41cl) . Consequently, the domain 
wall 


m{x) = m{y) := sgn{y), 
at y = 0, binds the zero mode 


''’^020 '’^030 1 ) I™ ,,i 

e 0 ■^ = e-\^^y\ 


G K, 


X, 


where 


^010 X = -sgn(m<,^)x 


(3.42a) 


(3.42b) 

(3.42c) 


with X independent of x and z. The kinetics of the gap¬ 
less boundary states is governed by the Dirac Hamilto¬ 
nian 

■^bd (2;, z) = -id^Tg O CTg - iS^Ti 0 CTg, (3.43) 

where we have chosen < 0. 

On the boundary, the operations for reversal of time 
and charge conjugation are now represented by 

:='^2 ® ^2 0 IK, :=iro 0(72 0 IK, 

(3.44a) 

where 1 is the v x v unit matrix. We seek the single¬ 
particle Hamiltonian on the boundary that encodes the 
fluctuations arising from the Hubbard-Stratonovich de¬ 
coupling of quartic interactions through a generic dy¬ 
namical mass that respects the PHS on the boundary. 
It is given by 

(D X. Z) ■■= - 0 (Tg 0 1 - id^T^ 0 (Tg 0 1 

-I- T 2 0 M(t, (c, z), (3.44b) 

with the 2v x 2v Hermitian matrix 

M{t, X, z) = +(a 2 0 1) M*{t, X, z) ((72 0 1). (3.44c) 

The space of normalized Dirac mass matrices satisfying 
the condition (I3.44cp is topologically equivalent to the 
space 


K := U ^Pi’^ - ^)] ( 3 - 45 ) 

fe=i 

for the symmetry class C in two-dimensional space. The 
limit cx) of these spaces is the classifying space R^. 
In order to gap out dynamically the boundary zero modes 
without breaking the defining symmetries of the symme¬ 
try class Cl, we need to construct a (2-|-l)-dimensional 
QNLSM for the (boundary) dynamical Dirac masses from 
the two-dimensional symmetry class C without topolog¬ 
ical obstructions. We construct explicitly the spaces for 
the relevant normalized boundary dynamical Dirac mass 
matrices of dimension v = 2'^ with n = 0,1,2 in the 
following. The relevant homotopy groups are given in 
Table 1^ 

Case V = \: There is one 2x2 Hermitian matrix 
M(t, X, z) on the boundary that is proportional to (7g. A 
domain wall in imaginary time such as m 2 00 sigii('r) T2 0 
(7 q prevents the dynamical generation of a spectral gap 
on the boundary. 

Case V = 2: We use the representation 1 = ctq. The 
Hermitian 4x4 matrix M(t,x,z) is a linear combina¬ 
tion with real-valued coefficients of the matrices X,,,,, = 
® cr/, with p, p' = 0,1,2, 3 such that A2g A20 = 
-l-A^^,. Of these, one finds the five matrices A^2; 
A22, A32 , Aq 3, and Agg that anticommute pairwise. If 
M(t,x,z) is a linear combinations with real-valued co¬ 
efficients of these five matrices and if M{t,x,z) is nor¬ 
malized by demanding that it squares to Agg, then the 
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TABLE VI. Reduction from Z to Z 4 for the topologically 
equivalent classes of the three-dimensional SPT phases in the 
symmetry class Cl that arises from interactions. We denote 
by V,^ the space oi v x v normalized Dirac mass matrices 
in boundary (d = 2) Dirac Hamiltonians belonging to the 
symmetry class C. The limit 1 / —>■ 00 of these spaces is the 
classifying space R^. The second column shows the stable D- 
th homotopy groups of the classifying space R 4 . The third 
column gives the number v of copies of boundary (Dirac) 
fermions for which a topological obstruction is permissible. 
The fourth column gives the type of topological obstruction 
that prevents the gapping of the boundary (Dirac) fermions. 


D 

WD{Ri) 

V 

Topological obstruction 

0 

Z 

1 

Domain wall 

1 

0 



2 

0 



3 

0 



4 

z 

2 

WZ term 

5 

Z2 

4 

None 

6 

Z2 



7 

0 




set spanned by M{t,x,z) is homeomorphic to 5^. As 
7r2+i+i(5'^) = Z, it is then possible to add a topological 
term to the (2-|-l)-dimensional QNLSM on the bound¬ 
ary that is of the WZ type. This term is conjectured to 
prevent the gapping of the boundary states. 

Case = 4.- We use the representation 1 = CTq ® ctq . 
The Hermitian 8x8 matrix M (t, x, z) is a linear com¬ 
bination with real-valued coefficients of the matrices 
cr^' ® ^’1," with /r" = 0,1, 2, 3 such 
that ^^ 200^200 = Of these, one finds 

the six matrices X 120 , ^ 220 ; ^3201 ^010 j ^031 ^033 

that anticommute pairwise. If M(t, x, z) is a linear com¬ 
binations with real-valued coefficients of these six matri¬ 
ces and if M{t,x,z) is normalized by demanding that 
it squares to ATqqq, then the set spanned by M{t,x,z) 
is homeomorphic to . It is then impossible to add a 
topological term to the (2-|-I)-dimensional QNLSM on 
the boundary. The boundary zero modes can be gapped 
out. 

We conclude that the effects of interactions on the 
three-dimensional SPT phases in the symmetry class Cl 
are to reduce the topological classification Z in the non¬ 
interacting limit down to Z4. The logic used to reach 
this conclusion is summarized by Table ED once the line 
corresponding to zz = 1 has been identified. It is given 
by the smallest D that accommodates a non-trivial entry 
for the corresponding homotopy group. 


3. The symmetry class AIII when d — 3 

By omitting the contributions arising from the gauge 
potentials, the single-particle Hamiltonian (I3.36al) does 
not specify uniquely the symmetry class. For example, 


the single-particle Hamiltonian (I3.36a|) can also be inter¬ 
preted as an insulator belonging to the symmetry class 
Am, for it anticommutes with the composition 

Ts := -iTC = A21 (3.46) 

of the operations T~ and C for time reversal and charge 
conjugation, respectively, defined in Eq. (j3.36d() . 

The direct product of the single-particle Hamiltonian 
(I3.36al) with the v x v unit matrix 1 supports u zero 
modes bound to the boundary ?/ = 0, for they are anni¬ 
hilated by the boundary Hamiltonian 

z) '■= Tg ® a - i5^ Ti 0 1 (3.47a) 

that anticommutes with 

r(“^=T2 0a. (3.47b) 

The fate of these zero modes in the presence of 
fermion-fermion interactions is investigated in two steps. 
First, we include the effects of interactions by perturbing 
the boundary Hamiltonian with all boundary dynami¬ 
cal Dirac masses from the symmetry class A. Second, we 
introduce a BdG (Nambu) grading to account for the 
interactions-driven superconducting fluctuations by per¬ 
turbing the boundary Hamiltonian TL^^^^x, z) ® Pq with 
all boundary dynamical Dirac masses that anticommute 
with rQ0]l0p4 K, i.e., with all boundary dynamical Dirac 
masses from the symmetry class D. In the hrst step, the 
boundary dynamical single-particle Hamiltonian is 

X, z) := (-id^ Tg -id^T^)®l 

+ T 2 ® M(t,x, z), (3.48a) 

with the V X V Hermitian matrix 

M{t^x^z) = {t^x^z). (3.48b) 

In the second step, the boundary dynamical single¬ 

particle Hamiltonian is 

'^hdt '(d X, z) := {-\d^ Tg - id^ Tg) 0 (po 0 1) 

-\- T 2 ® M{t,x,z), (3.49a) 

with the 2v x 2v Hermitian matrix 

M{t, X, z) = +{pi 0 1) M*{t, X, z) {pi 0 1). (3.49b) 

The space of boundary dynamical Dirac mass matri¬ 
ces of the form (|3.48bp that square to the unit matrix is 
homeomorphic to the classifying space Cq for the sym¬ 
metry class A in two-dimensional space and in the limit 
iz —>• cx). In order to gap out dynamically the bound¬ 
ary zero modes without breaking the defining symme¬ 
tries of the symmetry class AIII, we need to construct a 
(2-|-l)-dimensional QNLSM for the (boundary) dynam¬ 
ical Dirac masses from the two-dimensional symmetry 
class A without topological obstructions. When v = 1, 
a domain wall such as M (r, x, z) = sgn(r) prevents 

















16 


the gapping of the boundary zero modes. When v = 2, 
we choose the representation 1 = (Tq. The set spanned 
by M{t,x, z) = with the real-valued 

functions m^- (r, x, z) obeying the normalization condition 
a;, z) = 1 supports a monopole that binds 
zero modes in (2-|-l)-dimensional space and time, as 
7 r2(S'^) = Z. When = 4, we choose the representation 
a = Xqq where := cr^ ig) cr^, for /i, fi' = 0,1, 2,3. We 

may then write M (r, x, z) = (t, x, z) . 

Any X^^, other than the unit matrix Aqq belongs to a 
multiplet of five pairwise anticommuting matrices of the 
form Xjj^, ^ A'qq. Hence, we may always construct a 
set of normalized M{t,x,z) homeomorphic to S'^. Since 
7 r2_|_i_|_i = Z, it is possible to augment the corre¬ 

sponding boundary QNLSM in (2-|-l)-dimensional space 
and time by a WZ term that modifies the equations of 
motion in a local way. This term is conjectured to pre¬ 
vent the gapping of the boundary states. When v = 2'^ 
with n > 3, we choose the representation 1 = Aqq... 
where X^^,... ■= a^ ® ® ■ ■ ■ for = 0,1, 2, 3. 

Any X^^i,,, other than the unit matrix Aqq... belongs to 
a multiplet of no less than seven pairwise anticommut¬ 
ing matrices. It is for this reason that the boundary 
states are then necessarily gapped, for it is not permissi¬ 
ble to add a topological term to the action of the bound¬ 
ary QNLSM for a sphere of dimension larger than four. 
We conclude that the effects of interactions in the three- 
dimensional SPT phases in the symmetry class AIII is 
to reduce the topological classification Z in the noninter¬ 
acting limit down to Zg under the assumption that only 
fermion-number-conserving interacting channels are in¬ 
cluded in the stability analysis. The logic used to reach 
this conclusion is summarized by Table IVTIl once the line 
corresponding to v = 1 has been identified. It is given 
by the smallest D that accommodates a non-trivial entry 
for the corresponding homotopy group. 

The space of boundary dynamical matrices that sat¬ 
isfy the condition (|3.49bp and square to the unit ma¬ 
trix is homeomorphic to the classifying space i?2 for the 
symmetry class D in two-dimensional space and in the 
limit ly —>■ oo. Because the dimension of the boundary 
dynamical matrix (I3.49b[) is twice that of the boundary 
dynamical matrix (|3.48bl) . one might have guessed that 
the gapping of the boundary zero modes takes place for 
a value of v smaller than eight. This is not so, how¬ 
ever, because of two constraints. The first constraint is 
that of PHS. The second constraint restricts the target 
space for the boundary QNLSM that is built out of the 
boundary dynamical Dirac masses. The target space of 
the QNLSM must be invariant as a set under the action 
of a global gauge U{1) transformation that is generated 
by Tq ® P 3 ® 1. This global C/(l) symmetry implements 
conservation of the fermion number. Indeed, one verifies 
the following facts. When 1 ^ = 1, the boundary dynam¬ 
ical matrix M{t,x,z) is a linear combination of and 
P 2 with real-valued functions as coefficients. Hence, the 
space of normalized boundary dynamical Dirac mass ma- 


TABLE VII. Reduction from Z to Zg for the topologically 
equivalent classes of the three-dimensional SPT phases in the 
symmetry classes AIII that arises from the fermion-number- 
conserving interacting channels. We denote by W the space 
oi u X u normalized Dirac mass matrices in boundary {d — 2) 
Dirac Hamiltonians belonging to the symmetry class A. The 
limit —>• 00 of these spaces is the classifying spaces Cq. The 
second column shows the stable D-th homotopy groups of the 
classifying space Cq. The third column gives the number v of 
copies of boundary (Dirac) fermions for which a topological 
obstruction is permissible. The fourth column gives the type 
of topological obstruction that prevents the gapping of the 
boundary (Dirac) fermions. 


D 

'^d(C'o) 

V 

Topological obstruction 

0 

Z 

1 

Domain wall 

1 

0 



2 

z 

2 

Monopole 

3 

0 



4 

z 

4 

WZ term 

5 

0 



6 

z 

8 

None 

7 

0 




trices is homeomorphic to and invariant as a set under 
any global gauge C/(l) transformation when v = 1. Be¬ 
cause of 7r]^(iS'^) = Z, vortex lines bind zero modes in 
(2-|-l)-dimensional space and time. When iz = 2, we 
represent the unit 4x4 matrix by Pq <8 ) ctq and we ex¬ 
pand any 4x4 Hermitian matrix as a linear combination 
with real-valued functions as coefficients of the sixteen 
matrices X^^, = p^® cr^, with p, p' = 0,1, 2,3. The 
boundary dynamical matrix M{t, x, z) is a linear combi¬ 
nation with real-valued functions as coefficients of the ten 
matrices A'qq, Xq^, Algg, A^q, X-^, A^g, A2Q, X 21 , A2g, 
Ag2 that satisfy the constraint A^^, = -l-Agg A*^, A^g. 
Other than the unit matrix Agg, any one of these nine 
matrices belongs to a triplet of pairwise anticommuting 
matrices. However, not all such triplets are closed under 
the global U(l) transformation defined by multiplication 
from the left and right (conjugation) with Agg. How¬ 
ever, there exists a triplet that is closed under conjuga¬ 
tion by Agg. For example, each element from the triplet 
of pairwise anticommuting matrices Agj^, Agg, A32 is in¬ 
variant under conjugation with Agg. Moreover, no other 
matrix, satisfying the condition A^^, = -l-Agg A*^, Aj^g, 
anticommutes with this triplet. Hence, this triplet spans 
a set of normalized boundary dynamical Dirac masses 
that is homeomorphic to S^, each point of which is in¬ 
variant under the global C/(l) transformation associated 
with the conservation of the fermion number. Because 
of 7I'2(S'^) = Z, monopoles bind zero modes in (2-1-1)- 
dimensional space and time that prevent the gapping 
of the boundary states when v = 2. When ^ = 4, 
we represent the unit 8x8 matrix by Pg 0 ag ® dg 
and we expand any 8x8 matrix as a linear combina¬ 
tion with real-valued coefficients of the sixty-four matri- 
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ces = P/x ® V ® = 0,1,2,3. 

The boundary dynamical matrix M{t,x,z) is a linear 
combination with real-valued functions as coefficients of 
those matrices = +^ioo ^;!I/x'm" ^loo- ^ Other 

than the unit matrix ^oooi ot those matrices be¬ 

long to a quintuplet of pairwise anticommuting matrices. 
Among these, each element from the quintuplet ^ooii 
Aqqj, X^i 2 j Xq 22 i -^332 i® invariant under conjugation 
by AgQQ. Moreover, no other matrix, satisfying the con¬ 
dition = TAioo Aioo, anticommutes with 

this quintuplet. Hence, this quintuplet spans a set of nor¬ 
malized boundary dynamical Dirac masses that is homeo- 
morphic to S '^, each point of which is invariant under the 
global U{1) transformation associated with the conserva¬ 
tion of the fermion number. Because of 'n 2 +i+i{S'^) = Z, 
it is possible to augment the corresponding boundary 
QNLSM in (2-|-l)-dimensional space and time by a WZ 
term that modifies the equations of motion in a local 
way. This term is conjectured to prevent the gapping of 
the boundary states. When v = 2"“^ with n > 4, any 
permissible matrix X^^,^,„,„ = +X^qq... X*^^,X^^^... 
belongs to a (:/)-tuplet of pairwise anticommuting per¬ 
missible matrices with N{v) > 5i^ However, not all such 
A(:/)-tuplets are closed under the global [/(I) transfor¬ 
mation defined by conjugation with Aggp.... The N{v)- 
tuplet that contains the pair of anticommuting matri¬ 
ces AQg...gg = +Agg...gg AQg...g]^ -^10'"00 ^00'"03 = 

-|-AiQ...Qg Aqo...o 3 AiQ...Qg has the particularity that each 
of its elements is invariant under conjugation with Aggp... 
and cannot be augmented by one more anticommuting 

Hence, this N{^)- 

tuplet spans a set of normalized boundary dynamical 
Dirac masses that is homeomorphic to each 

point of which is invariant under the global C/(l) trans¬ 
formation associated with the conservation of the fermion 
number. Since N{i/) >5 for z/ = 2”“^ with n > 4, it fol¬ 
lows that all homotopy groups of order less than four for 
the space of the normalized boundary dynamical Dirac 
masses that are invariant under the global C/(l) trans¬ 
formation are vanishing. The boundary states are then 
necessarily gapped. 

We conclude that the effects of interactions on the 
three-dimensional SPT phases in the symmetry class AHI 
are to reduce the topological classification Z in the non¬ 
interacting limit down to Zg. 


4- The symmetry class All when d — 3 

We close the discussion of the stability to fermion- 
fermion interactions of strong noninteracting topological 
insulators or superconductors in three-dimensional space 
by illustrating why the Z2 noninteracting classification is 
stable. 

To this end, consider the single-particle bulk Dirac 


TABLE VIII. Stability to fermion-fermion interactions of 
the noninteracting topological classification Z2 for three- 
dimensional strong topological insulators belonging to the 
symmetry classes AIL We denote by W the space oi u x u 
normalized Dirac mass matrices in boundary (d = 2) Dirac 
Hamiltonians belonging to the symmetry class A. The limit 
zz —>■ 00 of these spaces is the classifying space Cq. The sec¬ 
ond column shows the stable D-th homotopy groups of the 
classifying space Cq. The third column gives the number u of 
copies of boundary (Dirac) fermions for which a topological 
obstruction is permissible. The fourth column gives the type 
of topological obstruction that prevents the gapping of the 
boundary (Dirac) fermions. 


D 

rrD(C'o) 

V 

Topological obstruction 

0 

Z 

1 

Domain wall 

1 

0 




Hamiltonian 

n^°\x) ■= -18^X21 -il 9 yAgi -ia,, Ag 2 -f TO(a;) Agg, 

(3.50a) 

where A^^^ := for /i, n' = 0,1, 2, 3. Because 

n‘-°\x)=+TV.^°\x)r-^, r:=iA2gK, (3.50b) 

we interpret this Hamiltonian as realizing a noninteract¬ 
ing topological insulator in the three-dimensional sym¬ 
metry class AH. The domain wall in the mass 

m{x, y, z) = TO^o sgn(2:) (3.51a) 

binds a zero mode to the boundary z = 0 that is annihi¬ 
lated by the boundary single-particle Hamiltonian 

^b°d ix,y) = - cr2 - idy CTg 

= (3.51b) 

where 

%d ■= K. (3.51c) 

The boundary dynamical Dirac Hamiltonian 

(d X, y) = -id^ 02 - \-dy o-^+M (r, x, y) erg (3.52) 

belongs to the symmetry class A, as the Dirac mass 
McTg breaks TRS unless M(—t, x, y) = —M{t,x, y). The 
space of normalized boundary dynamical Dirac mass ma¬ 
trices {icTg} is homeomorphic to the space of normalized 
Dirac mass matrices for the two-dimensional system in 
the symmetry class A 

1 

= U ^(^)/ [^(^) X - ^)] • (3.53) 

fe=0 

The domain wall in imaginary time M{T,x,y) = 
sgn(r) prevents the gapping of the boundary zero 
modes. 
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We conclude that the noninteracting topological clas¬ 
sification Z2 of three-dimensional insulators in the sym¬ 
metry class All is robust to the effects of interactions un¬ 
der the assumption that only fermion-number-conserving 
interacting channels are included in the stability analy¬ 
sis. The logic used to reach this conclusion is summa¬ 
rized by Table IVIIII once the line corresponding to v = 1 
has been identified. It is given by the smallest D that 
accommodates a non-trivial entry for the corresponding 
homotopy group. Moreover, one verifies by introducing 
a BdG (Nambu) grading that this robustness extends 
to interaction-driven dynamical superconducting fluctu¬ 
ations. 


D. Higher dimensions 

By working out explicitly the effects of fermion-fermion 
interactions on the boundary states supported by single¬ 
particle Dirac Hamiltonians representing strong topolog¬ 
ical insulators and superconductors when the dimension¬ 
ality of space ranges from d = 1 to d = 8, the following 
rules can be deducedi^ 

Rule 1: The Z2 topological classification of strong 
topological insulators and superconductors is robust to 
interactions in all dimensions. 

Rule 2: The Z topological classification of strong topo¬ 
logical insulators and superconductors is robust to inter¬ 
actions in all even dimensions. 

Rule 3: The Z topological classification of strong topo¬ 
logical insulators and superconductors is unstable to in¬ 
teractions in all odd dimensions. 

We prove Rules 1 and 2 in Secs. IIIID II and IIIID2I 
respectively. Finally, we work out explicitly the reduction 
pattern of the noninteracting Z topological classification 
for any odd dimension in Sec. IIIID 31 


1. The case 0/Z2 classification 

The proof of Rule 1 follows the same logic as in the 
example of the three-dimensional symmetry class All in 
Sec. lIIICT^ When d = 1, there are two symmetry classes 
with TT^iy) = Z2, the symmetry classes D and Dill (see 
Table |T|. No dynamical Dirac mass is allowed in class 
D, since there is no protecting symmetry to break. For 
the symmetry class Dill, the normalized boundary dy¬ 
namical Dirac masses are taken from Dirac masses in the 
symmetry class D and belong to the classifying space 
i?2: according to Table HI According to Table IXVil in Ap¬ 
pendix m 7ro(i?2) = ^2- Hence, the two noninteracting 
Z2 topological classification are stable in one-dimension. 
To treat the case of d > 2, let F denote any one of the 
eight real classifying spaces V and observe that, accord¬ 
ing to Table IXVil at least one of the homotopy groups 
ttd{V) with D = 0,1,2, 3 is non-trivial. We specialize 
to any one of the two symmetry classes in d dimensions 
with the classifying space (the space of normalized bulk 


Dirac masses) V such that 7rg(F) = Z2. By assumption 
d -I- 1 >3. Let I4,(j denote the space of the boundary 
dynamical Dirac masses. If this space is empty, the Z2 
classification is stable. If this space is not empty, then we 
know that at least one of with D = 0,1, • • • , d-bl 

is nonvanishing. In turn, this implies that at least one 
of the homotopy groups from Eq. (12.101) is non-trivial. 
As the sphere entering Eq. (12.101) is the target 

space for the QNLSM in (d — 1) -I- 1 space and time di¬ 
mensions obtained from integrating the v = 1 bound¬ 
ary Dirac fermions subjected to dynamical masses, the 
QNLSM accommodates a topological term that prevents 
the gapping of the v = \ boundary zero mode. Hence, 
the two noninteracting Z2 topological classification are 
stable in any spatial dimension. 


2. The case of even dimensions 

Because of Rule 1, we only need to consider the sym¬ 
metry classes in even dimensions which have Z topologi¬ 
cal classihcation for gapped noninteracting fermions. Ac¬ 
cording to Table [T] and the Bott periodicity of two (eight) 
for the complex (real) symmetry classes, these are the 
symmetry classes (i) A for d = 0 mod 2, (ii) AI and All 
for d = 4, 8 mod 8, (hi) D and C for d = 2,6 mod 8. 

Proof for case (i): We start with the complex sym¬ 
metry class A in even dimensions. We proceed in two 
steps. First we rule out dynamical superconducting fluc¬ 
tuations. We then show that the inclusion of dynamical 
superconducting fluctuations is harmless. 

Without dynamical superconducting fluctuations, the 
classifying space for the normalized dynamical Dirac 
masses is that for the complex symmetry class A. Because 
there is no symmetry that is violated by such dynamical 
Dirac masses, dynamical Dirac masses are forbidden al¬ 
together. 

With dynamical superconducting fluctuations, normal¬ 
ized dynamical Dirac masses form the space of Dirac 
masses in the symmetry class D. The original single¬ 
particle Hamiltonian that annihilates v zero modes 
is extended to a BdG (Nambu) single-particle Hamilto¬ 
nian that commutes with pg and anticommutes 

with pgK. Here, pg is the unit 2x2 matrix and p are 
the Pauli matrices acting on the auxiliary particle-hole 
degrees of freedom. Boundary dynamical Dirac masses 
may then exist. However, they must anticommute with 
P3, since no boundary Dirac mass that commutes with pg 
is allowed after restricting to the boundary. Upon 

integrating the boundary Dirac fermions, a QNLSM in 
(d — 1) -|-1 space and time dimensions ensues. The target 
space of this QNLSM has to be closed under the action of 
a global C/(l) gauge transformation. This is to say that 
a generic boundary dynamical Dirac mass must be of the 
form 


7' = [cos(0)pg -I- sin(0)p2] ® M, 


(3.54a) 
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with the J^/2) x (rjjjjjj z//2) matrix M satisfying 

M = = (3.54b) 

We recall that is the minimal rank of the BdG 
Hamiltonian. The action on 7' of a global U{1) gauge 
transformation parametrized by the global phase a is 
simply the shift 9 ^ 9 + a. If so, the target space of 
the QNLSM is homeomorphic to x ^ whereby 
is parametrized by 9 and is parametrized by M 

squaring to the unit matrix. For such a target space, we 
can always assign a topological term to account for the 
vortices supported by the parameter 9 for the factor, 
as 7r2(5'^) = Z. These vortices bind 1 / zero modes. 

We conclude that the noninteracting topological clas¬ 
sification Z for the symmetry class A in even dimensions 
survives strong interactions on the boundary provided 
the fermion-number conservation is neither explicitly nor 
spontaneously broken. 

Proof for case (ii): First, we show the statement for: 

(a) cases with dynamical Dirac masses that preserve the 
fermion-number U(l) symmetry. We then proceed to: 

(b) cases with C7(l)-breaking dynamical Dirac masses. 

(a) We consider the massive Dirac Hamiltonian 


algebra. We arrive at the Clifford algebra 

«3,d:={^7o,'r,Jr;7i,--- ,7d} (3-57) 

for the symmetry class AH. 

In both symmetry classes the choice of 7 q is unique, 
up to a sign, as a consequence of the fact that the zeroth 
homotopy groups of the classifying spaces for the sym¬ 
metry classes AI and AH is Z in 4n dimensions. In other 
words, no other Dirac mass matrix that is invariant un¬ 
der reversal of time anticommutes with 7 qi^ This leaves 
open the possibility that the Clifford algebras (13.561) and 
(13.571) for d = An could admit the addition of a gener¬ 
ator 7 q that anticommutes with 77^°^ and is odd under 
reversal of time, Jq = —T If so, the choice of 7^ 
to 7^ in the Clifford algebras p.56l) and (13.571) would not 
be unique in an uncountable (in a continuous) way. The 
existence of 7 q is thus tied to the task of parametrizing 
in a continuous way the representation of the generator 
(e.g., 7^) present in but absent in Clp ^ applied 

to the cases {p,q) = (l,4n-|-l) and {p,q) = (3,4n—I) for 
the 4n-dimensional symmetry classes AI and AH, respec¬ 
tively!^ Both tasks are denoted by the extension problem 
of Clifford algebras 


d 

TL^^^x) +m{x)'-fQ, m{x) gR, (3.55) 

1=1 


with the classifying spaces 


(3.58) 


obeying the TRS represented by T for classes AI and 
AH in d = 4n for n = 1,2, • • •. Here, the Dirac matrices 
are of dimension r > They obey the Clifford alge¬ 
bra {7^, 7^/} = 2 5^^, with pL, p' = 0, • • • ,d. The Dirac 
matrices entering %'A^\x) and the operator T that repre¬ 
sents reversal of time can be used to define the following 
pair of Clifford algebrasi^^*^ 

For the symmetry class AI, reversal of time is rep¬ 
resented by an element of the Clifford algebra T that 
squares to the unit matrix. It and iT, together with 
the gamma matrices 7i, ■ • ■, 7^ satisfying the conditions 
7i = -T7i T“^, • • •, 7(i = —T 7d'r”^, are generators of 
the Clifford algebra. On the other hand, the Dirac mass 
matrix iyg is chosen as the generator that squares to mi¬ 
nus the unit matrix, for 7g = +T 7gT“^. We arrive at 
the Clifford algebra 

«i,2+d:={^7o;'r,Jr,7i,--- ,7j (3.56) 

for the symmetry class AI, where J is the generator that 
represents the imaginary unit “i” and satisfies the rela¬ 
tions = — 1 and {T, J} = Oi^ 

For the symmetry class AH, reversal of time is rep¬ 
resented by an element of the Clifford algebra T that 
squares to minus the unit matrix. It and VT enter on 
equal footing with iyg as the generators that square to 
minus the unit matrix, for 7g = -|-7~7g T~^■ On the other 
hand, gamma matrices 7]^, • • •, 7^, satisfying the condi¬ 
tions 7i = •••, 7^ = are the 

generators that square to the unit matrix in the Clifford 


R 


q-p 


(p,?) = (l,4n +1), 
Rau-a, (p,?) = (3,4n-1), 


(3.59) 


as solutions for the set of representations of possible 7g in 
the symmetry classes AI and AH, respectively. Hereto, 
it is the zeroth homotopy group of the classifying space 
Rp_g that seals the fate of the existence of 7g. As 
7ro(i?4„) = 7rg(i?4„_4) = Z, it follows that 7g does not 
exist, i.e., no dynamical Dirac mass that breaks the TRS 
symmetry but preserves the global t/(l) gauge symmetry 
is permissible for the symmetry classes AI and AH when 
d = 4n. 

(b) After having established the absence of C/(I)- 
preserving dynamical Dirac masses, the stability anal¬ 
ysis in the presence of dynamical superconducting fluc¬ 
tuations for the symmetry classes AI and AH is the same 
as that for the symmetry class A. The boundary dynam¬ 
ical Dirac mass takes the form (13.541) . The target space 
of the QNLSM is homeomorphic to x since it 

has to be closed under the action of a global U{1) gauge 
transformation. Vortices that bind boundary zero modes 
originate from the manifold. 

Proof for case (Hi): The symmetry classes D and C for 
d = 2,6 (mod 8) do not support dynamical Dirac masses 
along the boundary, because the PHS is kept as a funda¬ 
mental symmetry. Their noninteracting topological clas¬ 
sification Z survives strong interactions on the boundary 
provided the PHS is neither explicitly nor spontaneously 
broken. 
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TABLE IX. Application of the Bott periodicity obeyed by the homotopy groups for D — 0,1, ■ ■ ■ of a given classifying 

space V^_i of dynamical boundary Dirac masses to deduce the reduction pattern Z —>■ for the symmetry class BDI in 

dimensions (a) d = 1 for which V'j_i = R 2 and i^max = 8, (b) d = 5 for which = Re and r'max = 16, and (c) d = 9 for which 
= i ?2 and i^niax — 128. The column v fixes the rank r := rjj,in i' of the Dirac Hamiltonian in the symmetry class BDI. The 
fourth column gives the target manifold of the QNLSM with the action S'qi.jlsm th^-f encodes the fermion-fermion interactions 
on the (d — l)-dimensional boundary. The fifth column indicates if a topological action Stop can be added to S'qj.^lsm- 

(a) Symmetry class BDI in d = 1 (b) Symmetry class BDI in d = 5 (c) Symmetry class BDI in d = 9 


D 710(8.2) 

u 

‘S'qnlsm 

‘ 5 'top 

D 

T^DiRe) 

p 

‘S'qnlsm 

‘ 5 'top 

D 

770(82) 

p 

S'qnlsm 

‘^top 

0 

Z2 

2 


/ 

0 

0 




0 

Z2 

2 

S° 

/ 

1 

0 




1 

0 




1 

0 




2 

z 

4 


/ 

2 

z 

1 


/ 

2 

z 

4 


/ 

3 

0 




3 

Z2 

2 


/ 

3 

0 




4 

0 




4 

Z2 

4 


/ 

4 

0 




5 

0 




5 

0 




5 

0 




6 

z 

8 


- 

6 

z 

8 


/ 

6 

z 

8 

s® 

/ 

7 

Z2 

16 


- 

7 

0 




7 

Z2 

16 

5 " 

/ 






8 

0 




8 

Z2 

32 


/ 






9 

0 




9 

0 









10 

z 

16 

^10 

- 

10 

z 

64 

510 

/ 











11 

0 














12 

0 














13 

0 














14 

z 

128 

514 

- 


3. The case of odd dimensions 

The topological classification Z of noninteracting 
strong topological insulators and superconductors in odd 
dimensions of space is reduced to the coarser classifica¬ 
tion Z^ with t'jnax integer: 


and Cl are interpreted as superconductors. 
Step f: Assign the minimal value 

'1, ^ora = 0, 

^min • ^ 

.2, 


(3.60b) 


. (3.60a) 

*^max ' ' 

The label “max” stands here for maximum. The task at 

hand is thus to compute the integer r'jjjg.x- Computing 

^max proceeds with the following algorithm (see Tables 

mm . 

Step 1: Choose any one of the ten AZ symmetry classes 
from Table ID 

Step 2: Choose any odd dimension d for which the zero¬ 
th homotopy group of the classifying space of the chosen 
symmetry class is the set of integers [7rQ(V)j) = Z]. 
This step restricts the symmetry classes to the complex 
symmetry class AIII and the real symmetry classes BDI, 
Dili, CII, and CL 

Step 3: Identify the parent symmetry class and its clas¬ 
sifying space Dj that follows if CHS is broken for the 
complex symmetry class AIII or if TRS is broken for 
the real symmetry classes. This step restricts the parent 
symmetry classes to the complex symmetry class A if the 
symmetry class AIII is interpreted as realizing an insu¬ 
lator, the real symmetry class D if the symmetry classes 
BDI and Dill are interpreted as superconductors, and 
the real symmetry class C if the symmetry classes CII 


if the zero-th homotopy group of is trivial or non¬ 
trivial, respectively. 

Step 5: Identify the classifying space that deter¬ 
mines the dynamical Dirac mass matrices induced by the 
fermion-fermion interactions on the boundary. 

Step 6: Construct a table with lines labeled by the in¬ 
teger D = 0,1,2,---. The first column gives the order 
D of the homotopy group given in the second 

column. The third column is the number v of bound¬ 
ary zero modes in the selected symmetry class. Enter 
the value third column for the smallest value 

of D for which 7r^(V)(_i) is non-trivial. The value of v 
is then doubled for each successive line with 7r^(V)(_i) 
non-trivial. The fourth column denotes the target space 
of the QNLSM with the action S'qnlsm defined by inte¬ 
grating out all the boundary Dirac fermions when cou¬ 
pled to (D -t- I) real-valued bosonic fields, each one of 
which couples to a Dirac mass matrix from a [D -|- 1)- 
tuplet of pairwise anticommuting Dirac mass matrices 
allowed on the boundary by the parent symmetry class. 
The fifth column indicates when a topological term 
can be added to the action 
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TABLE X. Application of the Bott periodicity obeyed by the homotopy groups iToiV) for -D = 0,1, • • • of a given classifying 
space Vd-\ of dynamical boundary Dirac masses to deduce the reduction pattern Z ^ for the symmetry class Dill in 

dimensions (a) d = 3 for which Vd-\ = Rq and i^max = 16, (b) d = 7 for which V^-i = R 4 and i/max = 32, and (c) d = 11 for 
which Vd-i = -Rq and z^max — 256. The column 1 / fixes the rank r := 1 / of the Dirac Hamiltonian in the symmetry class 

DHL The fourth column gives the target manifold of the QNLSM with the action Sg^Lgiyi that encodes the fermion-fermion 
interactions on the (d — 1)-dimensional boundary. The fifth column indicates if a topological action Stop can be added to 
'S'qnlsm- 

(a) Symmetry class Dill in d = 3 (b) Symmetry class Dill in d = 7 (c) Symmetry class Dill in d = 11 


D ttu{Ro) 

^ 'S'qnlsm 

•^top 

D 

T^oiRi) 

u 

'5'qnlsm 

a 

0 

D TVjy{Ro) 

i/ 

'S'qnlsm 

•^top 

0 

Z 

1 

/ 

0 

z 

1 

S° 

/ 

0 

Z 

1 

S° 

/ 

1 

Z 2 

2 

/ 

1 

0 




1 

Z 2 

2 

S' 

/ 

2 

Z 2 

4 

/ 

2 

0 




2 

Z 2 

4 

s^ 

/ 

3 

0 



3 

0 




3 

0 




4 

z 

00 

/ 

4 

z 

2 

s" 

/ 

4 

z 

8 

s" 

/ 

5 

0 



5 

Z 2 

4 

s® 

/ 

5 

0 




6 

0 



6 

Z 2 

8 

s« 

/ 

6 

0 




7 

0 



7 

0 




7 

0 




8 

z 

16 S'® 

- 

8 

z 

16 

s® 

/ 

8 

z 

16 

s® 

/ 





9 

0 




9 

Z 2 

32 

s® 

/ 





10 

0 




10 

Z 2 

64 

Sio 

/ 





11 

0 




11 

0 








12 

z 

32 

512 

- 

12 

z 

128 

Sl 2 

/ 


13 0 

14 0 

15 0 


16 Z 256 5^® 


Step 1: Let be the number of lines with Cl as superconductors, or 

non-trivial when D takes the values D = - ■ ■ ,d+l. 


It then follows that 

d = 8n+ 1 

d = 8n -I- 3 

d = 8n + 5 

d = 8n + 7 

BDI 

Z — y ^24/1+2 

- 

Z — y Z24n + 3 

- 

'^max = '^minX2"wz. (3.60c) (.jj 

Z —>• Z24n + 1 

- 

Z — y Z24n+4 

- 


For the complex symmetry class AlII in dimension d = 
2n-|-1 with n = 0,1, 2, • • •, the reduction pattern induced 
by the fermion-fermion interactions is 

Z-5>Z2„+2. (3.61) 

By making use of the eight-fold Bott periodicity, one ver¬ 
ifies that the reduction patterns are 


(3.63) 

if we interpret the symmetry classes BDI and CII as in¬ 
sulators. We have thus found two different patterns for 
the reduction of the topological classification of the sym¬ 
metry class BDI depending on these two interpretations, 
as we have observed for d = 1 in Sec. IIII A 31 



d = Sn -I- 1 

d = 8n -I- 3 

d = 871 -I- 5 

d = 8n -I- 7 

BDI 

Z — y Z2471+3 

- 

Z — y Z24T1+4 

- 

Dili 

- 

Z — y Z24n+4 

- 

Z ^ Z24nd-5 

CII 

Z — y Z24n + 1 

- 

Z —y Z24n-f4 

- 

Cl 

- 

Z —y Z24n+2 

- 

Z — y Z24n + 5 


(3.62) 

if we interpret the symmetry classes BDI, Dill, CII, and 


IV. REDUCTION FOR TCS AND TCI 

The addition of discrete symmetries (such as spatial 
ones) to the PHS, TRS, and CHS enriches the classifica¬ 
tion of noninteracting fermions. A general method to ac¬ 
count for additional symmetries that square to the unity 
has been proposed in Ref. [l^. Hereto, local fermion- 
fermion interactions can reduce the noninteracting topo- 
fogical classification over Z, as was first observed in Refs. 
[Tr [13, Izlj and [tI by way of explicit examples. We are 
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TABLE XL Application of the Bott periodicity obeyed by the homotopy groups for D — 0,1, ■ ■ ■ of a given classifying 

space Vd_i of dynamical boundary Dirac masses to deduce the reduction pattern Z —>■ for the symmetry class CII in 

dimensions (a) d = 1 for which V'j_i = Rg and z^max = 2, (b) d = 5 for which = R 2 and t'max = 16, and (c) d = 9 for which 
= Rg and z^max = 32. The column 1 / fixes the rank r := rjj,in ^ of the Dirac Hamiltonian in the symmetry class CII. The 
fourth column gives the target manifold of the QNLSM with the action S'qj.^lsm that encodes the fermion-fermion interactions 
on the (d — l)-dimensional boundary. The fifth column indicates if a topological action Stop can be added to 

(a) Symmetry class CII in d = 1 (b) Symmetry class CII in d = 5 (c) Symmetry class CII in d = 9 


D TVniRa) 

U 

S'qnlsm 

‘^top 

D 7r£)(7?2) 

U 

'Sqnlsm 

‘^'top 

D 

T^DiRe) 

p 

‘Sqnlsm 

Stop 

0 

0 




0 

^2 

2 

5° 

/ 

0 

0 




1 

0 




1 

0 




1 

0 




2 

z 

1 


/ 

2 

z 

4 

S^ 

/ 

2 

z 

1 

S^ 

/ 

3 

^2 

2 

S^ 

- 

3 

0 




3 

Z 2 

2 

S^ 

/ 

4 

Z 2 

4 

S^ 

- 

4 

0 




4 

Z 2 

4 

S^ 

/ 

5 

0 




5 

0 




5 

0 




6 

z 

8 

5® 

- 

6 

z 

8 

5® 

/ 

6 

z 

8 

5® 

/ 

7 

0 




7 

Z 2 

16 

5" 

- 

7 

0 





8 0 

9 0 

10 Z 16 S^° / 

11 Z 2 32 


TABLE XII. Application of the Bott periodicity obeyed by the homotopy groups for Z) = 0,1, • • • of a given classifying 

space of dynamical boundary Dirac masses to deduce the reduction pattern Z ^ Z^ for the symmetry class Cl in 

dimensions (a) d = 3 for which = R^^ and z^max = 4, (b) d — 7 for which V^_i — Rg and z^max = 32, and (c) d = 11 for which 
V^_i — R 4 and t'max = 64. The column u fixes the rank r ~ v of the Dirac Hamiltonian in the symmetry class CL The 
fourth column gives the target manifold of the QNLSM with the action S'qnlsm that encodes the fermion-fermion interactions 
on the (d — l)-dimensional boundary. The fifth column indicates if a topological action Stop can be added to S'qnlsm- 

(a) Symmetry class Cl in d = 3 (b) Symmetry class Cl in d = 7 (c) Symmetry class Cl in d = 11 


D TVniRi) 

p 

'S'oNLSM 

•^top 

D TVniRo) 

p 

■SoNrsM 

0 

V 

D 

t^d{Ra) 

p 

'S'oNLSM 

•^top 

0 

Z 

1 


/ 

0 

Z 

1 

5® 

/ 

0 

Z 

1 


/ 

1 

0 




1 

Z2 

2 


/ 

1 

0 




2 

0 




2 

Z2 

4 


/ 

2 

0 




3 

0 




3 

0 




3 

0 




4 

z 

2 


/ 

4 

z 

8 


/ 

4 

z 

2 

5" 

/ 

5 

Z2 

4 

S® 

- 

5 

0 




5 

Z2 

4 

S'® 

/ 

6 

Z2 

8 

S® 

- 

6 

0 




6 

Z2 

8 

S® 

/ 

7 

0 




7 

0 




7 

0 









8 

z 

16 


/ 

8 

z 

16 

s« 

/ 






9 

Z2 

32 

5® 

- 

9 

0 














10 

0 














11 

0 














12 

z 

32 

^12 

/ 











13 

Z2 

64 

s« 

- 


interested in the robustness to local fermion-fermion in¬ 
teractions of noninteracting topological phases with re¬ 
flection symmetry. We treat two-dimensional topological 
superconductors in the symmetry class Dill with an ad¬ 
ditional reflection symmetry (the Yao-Ryu model from 
Ref. [H) and three-dimensional topological insulators in 


the symmetry class All with an additional reflection sym¬ 
metry (a topological crystalline insulator is realized in 
SnTe, as was shown in Ref. [zi- We show that the re¬ 
duction Z ^ Zg holds in both cases by applying the 
approach detailed in Sec. m 

The notation X , „ := t,, ® t', 0 t"„ 0 ••• for 

fXfX fx • • • /X /X /X 
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/i, /i', /i", • • • = 0,1, 2, 3 where Tq, Tq, Tq', • ■ • are unit 
2x2 matrices and the other r^, t"„, ••• are the 

corresponding Pauli matrices. 

A. Two-dimensional superconductors with 
time-reversal and reflection symmetries (Dlll+ii) 

Let X = {x,y) = {xi,X 2 ) denote a point in two- 
dimensional space. Let := (8> with = 

0,1, 2, 3 denote a basis for the vector space of all Hermi- 
tian 4x4 matrices. Following Yao and Ryuji^ consider 
the two-dimensional bulk single-particle Dirac Hamilto¬ 
nian 

'H^°\x,y) := -iS^Ygi - idyX^^ + m{x,y)XQ 3 . (4.1a) 

This single-particle Dirac Hamiltonian belongs to the 
symmetry class Dili, for 

rn^°\x,y)r-^ = +n^°Hx,y), (4.1b) 

Cn^°Hx,y)C-^ = -n^°\x,y), (4.1c) 

where 

T:=iX2oK, C:=XoiK. (4.1d) 

In addition, the Dirac Hamiltonian is invariant under re¬ 
flection in the x direction, 

7^, (-X, y ){ n ,)-^ = (x, y), (4.1e) 

where 

7^, = iX^o. (4. If) 

The operators Y, C, and TZ^ commute pairwise and 
square to 7^ = —1, = -1-1, and 72.^ = —1. 

The Dirac mass matrix Xq^ is here the only one al¬ 
lowed for dimension r = = 4 Dirac matrices under 

the symmetry constraints (I4.1bl) . (l4.1cL and (14. lei) . The 
domain wall 

m{x, y) = sgn(y) (4.2a) 


The boundary single-particle Hamiltonian (j4.2dl) and 
the operators (l4.2eD are denoted by 

and respectively, when tensored with the v x v 

unit matrix 1. The single-particle Hamiltonian 'H^)^^^{x) 
supports v linearly independent zero modes. Their sta¬ 
bility to interactions that preserve the symmetries is 
probed by studying the dynamical single-particle bound¬ 
ary Hamiltonian 

(d x) := O 1 -h -/'{t, x), (4.3a) 

where the boundary dynamical Dirac mass matrix {x) 
satisfies the particle-hole symmetry 

= -y (4.3b) 

and takes the form 

7'(x) = (72 ® x) + 7W2(t, x) (4.3c) 

with the v X v Hermitian matrices 

(d x) = +Ml (r, x ), M2 (r, x) = -M^ (r, x ), 

. 

i.e., Mi{t,x) is a real-valued symmetric matrix while 
M2 (r, x) is an imaginary-valued antisymmetric matrix. 
This is to say that the normalized boundary dynamical 
Dirac mass matrix 7'(a;) belongs to the space 

14 := 0{v), = lim 0{u), (4.3e) 

y—¥oo 

for the Dirac matrices in the boundary (d = 1) Dirac 
Hamiltonians belonging to the symmetry class D>2^ Inte¬ 
grating the boundary Dirac fermions delivers a QNLSM 
in (l-l-l)-dimensional space and time. In order to gap out 
dynamically the boundary zero modes without breaking 
the time-reversal, particle-hole, and reflection symme¬ 
tries, this QNLSM must be free of topological obstruc¬ 
tions. We construct explicitly the spaces for the relevant 
normalized boundary dynamical Dirac mass matrices of 
dimension ^ = 2” with n = 0,1, 2,3 in the following. The 
relevant homotopy groups are given in Table IXHIl 

Case ix = 1: There is a topological obstruction of the 
domain wall type as the target space is 


at y = 0 binds the zero mode 


= {±<Xy} (4.4a) 


-i^02 ^03 /dy 
e 0 = 


g-|™ooyl 


X, 


where 


^01 X = -sgn(moo) X 


(4.2b) 


(4.2c) 


with X independent of x and z, that is annihilated by the 
boundary Hamiltonian 

^b°d(a^) :=(4.2d) 

where we have chosen moo < 0. On the boundary y = 0, 
the symmetries (I4.1el) are realized by 

TOd := 1^2 K, Cbd := K, := i(T2. (4.2e) 


and 7rg(S'°) 7^ 0. 

Case V = 2: There is a topological obstruction of the 
vortex type as the target space is 

S'^ = {C1X21 -h C2X23I c^ + c^ = I, Ci e M} (4.4b) 
and 77^(5'^) = Z. 

Case 1 / = 4: There is a topological obstruction of the 
WZ type as the target space is 

='|c 3 X 23 o + C2X23Q -I- C3X3Q2 -|- C4X222 

^c-=I, c^gkI (4.4c) 

i=l 
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TABLE XIII. Reduction from Z to hg due to interactions 
for the topologically equivalent classes of the two-dimensional 
topological snpercondnctors protected by time-reversal and 
reflection symmetries (DIII-|-i?). We denote by the space 
oi V X V normalized Dirac mass matrices in boundary {d = 1) 
Dirac Hamiltonians belonging to the symmetry class D. The 
limit iz —> oo of these spaces is the classifying space Ri- The 
second colnmn shows the stable D-th homotopy groups of the 
classifying space . The third column gives the number v of 
copies of boundary (Dirac) fermions for which a topological 
obstruction is permissible. The fourth column gives the type 
of topological obstruction that prevents the gapping of the 
boundary (Dirac) fermions. 


D 

■kd{Ri) 

V 

Topological obstruction 

0 


1 

Domain wall 

1 

^2 

2 

Vortex 

2 

0 



3 

z 

4 

WZ term 

4 

0 



5 

0 



6 

0 



7 

z 

8 

None 


and = Z. 

Case V = 8.' There is no topological obstruction as one 
can find more than four pairwise anticommuting matrices 
such as the set 

{"^21001 ^2310) ^23311 ^2333) ^112o}- (4.4d) 

We conclude that the effects of interactions on the two- 
dimensional topological superconductors in the symme¬ 
try class Dill with additional reflection symmetry are to 
reduce the topological classification Z in the noninteract¬ 
ing limit down to Zg. 

B. Three-dimensional insulators with time-reversal 
and reflection symmetries (AII-I-R) 

We consider again the bulk, boundary, and dynamical 
boundary Hamiltonian defined in Sec. IIII C 41 i.e., Eqs. 
(13.501) (13.531) . We observe that the single-particle Hamil¬ 
tonian (13.50a[) has the symmetry 

(4.5a) 

where 

7^, :=iAio, 7^2 = -l, [r,7ej=0, (4.5b) 

in addition to the TRS (j3.50b(l . The presence of the 
additional reflection symmetry allows one to define a 
mirror Chern number (n_|_ € Z) for the sector with the 
eigenvalue = -(-i on the two-dimensional mirror plane 
(fcj, = 0) in the three-dimensional Brillouin zonei^ (The 
mirror Chern number for the eigensector TZ^ = —i is 


—n_|_.) Thus, the v linearly independent zero modes 
that follow from tensoring the single-particle Hamilto¬ 
nian (13.50all with the v x v unit matrix 1 along the do¬ 
main wall (13.51all are stable to strong one-body perturba¬ 
tions on the boundary that preserve the reflection sym¬ 
metry (14.51) . (If we forget the reflection symmetry and 
keep only the TRS, as we did in Sec. IHI C 41 it is only 
the parity of v that is stable to strong one-body pertur¬ 
bations on the boundary.) If we only consider dynamical 
masses that preserve the fermion-number U{1) symme¬ 
try, the space of normalized boundary dynamical Dirac 
mass matrices after tensoring the boundary dynamical 
Dirac Hamiltonian (13.531) with 1 is homeomorphic to the 
space of normalized Dirac masses for the two-dimensional 
system in the symmetry class A 

V,= \JU{,^)/[U{k)xU{:.-k)]. (4.6) 

/i-O 

The limit i/ ^ oo of these spaces is the classifying space 

Integrating the boundary Dirac fermions delivers a 
QNLSM in (2-|-I)-dimensional space and time. In order 
to gap out dynamically the boundary zero modes with¬ 
out breaking the symmetries, this QNLSM must be free 
of topological obstructions. We construct explicitly the 
spaces for the relevant normalized boundary dynamical 
Dirac mass matrices [M{T,x,y) in Eq. (13.521) ] of dimen¬ 
sion 1/ = 2" with n = 0,1, 2,3 in the following. The 
relevant homotopy groups are given in Table IXIVl 

Case ix = 1: There is a topological obstruction of the 
domain wall type as the target space is 

S° = {±1} (4.7a) 

and ^ 0. 

Case V = 2: There is a topological obstruction of the 
monopole type as the target space is 

5^ = {ciAi + C2A2 + C3A3I cl + cl + cl = 1} (4.7b) 

and 7r2(5'^) = Z. 

Case V = A: There is a topological obstruction of the 
WZ type as the target space is 

='|c3A33 -I- C2A23 -I- C3A33 -|- €4X33 -I- C5AQ2 

= (4.7c) 

and 7r4(5'"‘) = Z. 

Case 1 / = 8: There is no topological obstruction as one 
can find more than five pairwise anticommuting matrices 
such as the set 

{^133: ^ 233 j ^3331 ^013! ^023i ^OOlJ ^002}- (4-7d) 

We conclude that the effects of interactions on the 

three-dimensional topological insulators in the symme¬ 
try class AH with an additional symmetry are to reduce 
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TABLE XIV. Reduction from Z to Zg due to interactions for 
the topologically equivalent classes of the three-dimensional 
topological insulators with time-reversal and reflection sym¬ 
metries (AII-|-ii). We denote by the space of v xv normal¬ 
ized Dirac mass matrices in boundary {d = 2) Dirac Hamil¬ 
tonians belonging to the symmetry class A. The limit i/ —>■ oo 
of these spaces is the classifying space Cq. The second col¬ 
umn shows the stable D-th homotopy groups of the classifying 
space Cq. The third column gives the number v of copies of 
boundary (Dirac) fermions for which a topological obstruction 
is permissible. The fourth column gives the type of topolog¬ 
ical obstruction that prevents the gapping of the boundary 
(Dirac) fermions. 


D 

’’^oiCo) 

p 

Topological obstruction 

0 

Z 

1 

Domain wall 

1 

0 



2 

z 

2 

Monopole 

3 

0 



4 

z 

4 

WZ term 

5 

0 



6 

z 

8 

None 

7 

0 




the topological classification Z in the noninteracting limit 
down to Zg. 

This Zg classification is unchanged if all boundary 
dynamical masses that break the fermion-number i7(l) 
symmetry are accounted for. The corresponding tar¬ 
get spaces for the boundary dynamical masses and their 
topological obstructions are derived as was done in the 
stability analysis made for the symmetry class AIII in 
d = 3 in Sec. ImCSl Namely, we extend the single¬ 
particle Hamiltonian [Eq. (13.521) ] to a BdG Hamiltonian 

cr2 ® P3 - (g) Pg) (g) 1 -h (t, X, y), 

(4.8) 

where pg and p^ are unit 2x2 and Pauli matrices, re¬ 
spectively, acting on the particle-hole (Nambu) space and 
the particle-hole symmetry is given by C = Pi K. In 
this case, the target spaces of the QNLSM made of nor¬ 
malized boundary dynamical Dirac mass matrices 7' of 
dimension zz = 2" with n = 0,1,2,3 are modified as 
listed in the following with the notation = 

a^ (g> p^, (g) T^n (g) r^,„ . ... The relevant homotopy groups 
are given in Table DCVl We note that these target spaces 
are closed under the global 17(1) transformation gener¬ 
ated by P3. 

Case V = 1: There is a topological obstruction of the 
vortex type as the target space is 

= {c,X 2 i + C 2 X 22 I cf + cl = 1} (4.9a) 

and 7r3(5'^) = Z. 

Case 1 / = 2: There is a topological obstruction of the 


monopole type as the target space is 

S = {C1X21Q + C2A'22g + 03X332! Cl + C2 + = 1 } 

(4.9b) 

and 7r2(5'^) = Z. 

Case 1/ = 4.- There is a topological obstruction of the 
WZ type as the target space is 

=■!03X2333 -I- C2X2233 + C3X3323 -I- C4X3332 + C5X3332 
^0^ = 1, (4.9c) 

i=l ' 

and 7r4(5'^) = Z. 

Case v = 8: There is no topological obstruction as one 
can find more than five pairwise anticommuting matrices 
such as the set 

{^21000: ^22000) ^30200) ^30120) ^303121 ^30332}- (4.9d) 

Therefore, the topological classification Zg for three- 
dimensional TCIs in the symmetry class AH with an 
additional reflection symmetry is unchanged when the 
superconducting fluctuations are accounted for. This Zg 
classification is consistent with the results obtained re¬ 
cently in Refs, [t^ and [s^- We note that the classifying 
space Rq for the dynamical masses is the same as that in 
the case of three-dimensional TSs in the symmetry class 
DHL There is an important difference, however. Namely, 
the line corresponding to iz = 1 is moved to D = 1 in 
Table IXVI from D = 0 in Table |Vl This change origi¬ 
nates from the fact that the minimum matrix dimension 
of the BdG Hamiltonian for three-dimensional 

TGIs in the symmetry class All-l-i? is four, while that 
for three-dimensional TSs in the symmetry class Dili is 
two. Hence, the breakdown of the topological classifica¬ 
tion Z for three-dimensional TCIs in the symmetry class 
AH-fi? takes place at iz = 8, which is the half of zz = 16 
for three-dimensional TSs in the symmetry class DHL 


C. Massless Dirac fermions on the surfaces of SnTe 

The crystal SnTe is a three-dimensional topological 
crystalline insulator protected by time-reversal and re¬ 
flection symmetries (AH-fi?). SnTe supports four Dirac 
cones on the [001] surface and six Dirac cones on the 
[111] surface. If strong interaction effects are present, we 
expect that the zz = 4 phase described by a QNLSM with 
a WZ term should be realized on the [001] surface. At 
the [111] surface, we expect that the zz = 6 = 4 -f 2 phase 
be realized, whereby the effective field theory is that of 
a QNLSM with a WZ term for 4 out of the six surface 
Dirac cones and that of a QNLSM with a topological 
term arising from a gas of monopoles of the remaining 
two surface Dirac cones. 
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TABLE XV. Reduction from Z to Zg due to interactions for 
the topologically eqnivalent classes of the three-dimensional 
topological insulators with time-reversal and reflection sym¬ 
metries (All-fi?) when the superconducting fluctnations are 
accounted for. We denote by W the space of v xu normalized 
Dirac mass matrices in boundary {d — 2) Dirac Hamiltonians 
belonging to the symmetry class D. The limit —>■ oo of these 
spaces is the classifying space Rq. The second column shows 
the stable D-th homotopy groups of the classifying space Rq. 
The third column gives the number v of copies of boundary 
(Dirac) fermions for which a topological obstruction is per¬ 
missible. The fourth column gives the type of topological ob¬ 
struction that prevents the gapping of the boundary (Dirac) 
fermions. 


D 

^d{Ro) 

u 

Topological obstruction 

0 

Z 



1 

Z 2 

1 

Vortex 

2 

Z 2 

2 

Monopole 

3 

0 



4 

z 

4 

WZ term 

5 

0 



6 

0 



7 

0 



8 

z 

8 

None 
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Appendix A: Defining symmetries of strong 
topological insulators (superconductors) 

Dehne the many-body quadratic form 

H = J d‘^x J 

ij 

(Ala) 

where 

'H,j{x,y) = 'H*^{y,x), (Alb) 

and 

{iP^{t,x),^l{t,y)} = 6^^6{x-y) (Ale) 

are the only non-vanishing equal-time anticommutators. 

Time-reversal symmetry. Let K denote complex con¬ 
jugation. Define the time-reversal transformation by the 
antiunitary transformation 

T:=rK (A2a) 

that reverses time but leaves space unchanged by de¬ 

manding that 

(A2b) 


and 

t^^.(t,y)f-i (A2c) 

j' 

One verifies that 

fHf-^=H (A3a) 

if and only if 

Y.%,.n:^{x,y)T^} = n,^,{x,y). (A3b) 

ij 

Particle-hole (charge-conjugation) symmetry. Assume 
that 

Y,n,,{x,y)=Q. (A4) 

i 

Define the particle-hole transformation by the unitary 
transformation 

C := C (A5a) 

that reverses the sign of the fermion number 

^ = 0) := {x) -ip^ix) - I S{x = 0) (A5b) 

measured relative to the background of the fermion den¬ 
sity 1/2 but leaves space unchanged by demanding that 

C-i = C1' (A5c) 

and 

C tpj {t, y)C-^ -p], {t, y). (A5d) 

j' 

One verifies that 

CHC-^=H (A6a) 

if and only if 

®)- (^6b) 

Chiral symmetry. Assume that 

= (A7) 

i 

Define the chiral transformation by the antiunitary trans¬ 
formation 

S := 5 K (A8a) 

that reverses the sign of the fermion number 

n^{x) -^S{x = 0) := p)} (a;) ip^^x) -^S{x = 0) (A8b) 


f-i ^ ft 
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measured relative to the background of the fermion den¬ 
sity 1/2 but leaves space unchanged by demanding that 

5 -^ = 5 !' (A8c) 

and 

(A8d) 

j' 

One verifies that 

SHS-^=H (A9a) 

if and only if 

b 

The unitary symmetry under C is called charge conju¬ 
gation symmetry or particle-hole symmetry (PHS). The 
antiunitary symmetry under S is called the chiral symme¬ 
try (CHS). The antiunitary symmetry under T is called 
time-reversal symmetry (TRS). 


Appendix B: Tenfold way and classifying spaces 

In this appendix, we summarize the classification of 
gapped phases of noninteracting fermions in terms of the 


tenfold way. We also define the classifying spaces of nor¬ 
malized Dirac masses. The ten Altland-Zirnbauer (AZ) 
symmetry classes for Hermitian matrices are shown in 
Table H There, two complex and eight real symmetry 
classes are characterized by the presence or the absence of 
time-reversal symmetry (T), particle-hole symmetry (C), 
and chiral symmetry (T). Their presence is indicated 
by the sign entering the squared operators, = ±1 or 
= ±1, and by 1 for T. Their absence is indicated 
by 0. For each symmetry class and for any dimension 
d = 0,1, 2,... of space, the classifying space V/, which is 
the space of normalized Dirac masses, is given in the last 
column by labels to symmetric spaces. We list the ten 
relevant symmetric spaces and their homotopy groups in 
the stable homotopy regime in Table IXVIl The number 
N is related to the dimension r = of the Dirac 

matrices, i.e., A = 1,2,... is the number of copies of the 
minimal massive Dirac Hamiltonian of rank The 

stable homotopy regime refers to the limit N ^ 00 . Ac¬ 
cording to the Bott periodicity, the complex classifying 
spaces obey the periodicity condition 

T^D{Cq)=TTjj^2iCg), (g = 0,l), (Bl) 

and the real classifying spaces obey the periodicity con¬ 
dition 

T^oiRq) = T^D+siRq)^ (9 = 0 , ...,7). (B2) 
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TABLE XVI. Complex and real classifying spaces and their stable homotopy groups. Homotopy groups 7r^(y) for complex 
and real classifying spaces are periodic in D with periods of 2 and 8, respectively. 
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Classifying space V 

MV) MV) MV) MV) MV) MV) MV) MV) 

Co 

U^i^o{UiN)/[Uin) X UiN - n)]} 
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Rq 
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The symbol Tr represents tracing over the single-particle 
Hilbert space of the Dirac Hamiltonian with the Dirac ma¬ 
trices q: and /3 of dimension r. The saddle point equations 
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We denote with 11^-1 the area of the unit sphere , with 
k the length of the vector (lo, k), and with A the ultraviolet 
cutoff in {uj,k) space. The saddle-point equations reduce 
to the equation 
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It has the solution 


|<^/3| = i</>o(Ar), <j>o{Xr)>0. 
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